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ABSTRACT 


In Part I the primary longitudinal stability derivatives of a 
"eonventional"” airplane and a "typical" large, high speed, single 
rotor helicopter are simply compared, and fundamental similarities 
and differences discussed. Physical interpretation and meaning is 
given to various aspects of the character of the two types of air- 


Chae nia ie 


In Part II the influence of rotor flapping (tip path plane 
orientation and coning) on the stability derivatives of the whole 
helicopter is investigated, and then various stabilization schemes 
in which the flapping of the rotor is the basis for the feedback 
signal are explored. In particular are the effects of Oehmichen 
and Delta Three pitch-flap couplings on the static stability of the 
helicopter. It is found that of the techniques investigated and 
within the limits of validity of this work, only negative Delta 
Three (positive pitch-flap coupling) is effective in improving the 


longitudinal gust response characteristics in high speed flight. 
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INTRODUC TION 


This thesis iS in two parts. The first is a rudimentary cco@parison 
of the longitudinal gust response characteristics of a "conventional" air- 
plane with a "typical" large, high speed, single rotor helicopter. The 
second part consists of the investigation of various forms of helicopter 
stabilization techniques, especially those using rotor flapping as the 
sensed element. These two subjects, while related, are themselves topics 
of considerable reach, and their being treated together in a master's 
degree thesis is indicative of the depth they are explored. There is a 


reasoning, however, behind this overview type of approach. 


Helicopter pilots have always complained about the so called flying 
qualities of their craft in gusty weather. "They bloy around too much.” 
In order to find out why they do so, and in what way, it is expedient to 
relate the problem to one whose characteristics are better know - the 
airplane. The problem is commonly formulated as a set of differential 
equations with constant coefficients, and the characteristics of the 
problem are the signs and magnitudes of the coefficients - the stability 
derivatives. Trends, order of magnitudes, and signs of these derivatives 
indicate fundamental differences and similarities between the two types 
of aircraft and lead to insight as to why the helicopter responds more 
severely to gusts. The rudimentary form of the analysis is such that 


only these types of conclusion can, and are, drawn. 


The first part is thus a review of the derivation of (what are con- 
sidered to be) the important stability derivatives for aircraft gust 
response, and the reasoning behind several approximations to the full 
equations of motion with the author's two cents worth of physical 
interpretation and understanding added. It is, so to speak, the state- 


ment of the problem to be solved in the second part. 





In the second part, again, there is a good deal of review, the 
emphasis being on the physical interpretation and consequences of the 
mineines., Ihe 2nfiluence of rotor flapping on the stability derivatives 
of the whole helicopter, and hence the equations that govern how the 
pilot will be bounced around, is investigated. The impetus for this is 
to try to discover, in a simple, fundamental manner, whetner or not 
Suse Migaonon LGechniques involving the reduction of rotor flapping 
as a means to improve helicopter flying qualities are as "good" as they 
might seem from a more cursory glance at the problem. This is done by 
first deleting the effects of flapping in the expressions for the 
stability derivatives, and then by investigating feedback techniques 
such as Oehmichen and Delta Three pitch-flap Soup ines: Much time is 
Beent with these two schemes because, as is shown, they form the 
basis for several other stabilizer techniques which have been proposed. 
Finally, all of these results are contrasted with a conventional body 


mounted rate and attitude SAS (stability augmentation system). 


All of the analytical analyses are complimented by dynamic analog 
Computer responses. The data which is used to generate the helicopter 
models is presented in Appendix I, and serves to characterize the large, 
high speed, single rotor helicopter. Additionally, in the appendices are 
given a treatment of axis system conversion,and the definition of gust 


alleviation as it is used in this work. 





DESCRIPTION OF ANALYSIS 


The equations of motion used as a model of the aircraft in this 
work are the standard, small perturbation equations commonly used in 
stability analyses. They are extremely simplified by the approximations 
and assumptions that form the basis for this study - and result in a form 
of equations in which only general trends and important characteristics 
Gan be discovered. The justification for all of the simplification is 
eMeanveres, in finding order of magnitude trends that are characteristic 
of fixed or rotating wing aircraft, rather than the details that may 
vary considerably among aircraft of the same type. Some of the assumptions 
are purely for simplicity and others are assumptions that have been shown 
beebe Valid fOr an analysis or this sort, provided certain restrretions 
MaemtepGy in mind. in particular are the effects of high gain feedback 
fectons involving the rotor in the feedback loop. In such cases knew 
stability boundaries are avoided. Throughout the analysis the simplifying 
assumptions almost invariably are optimistic with regards to stability. 
Por this reason it 1S considered desirable to obtain results of a simple 


analysis before a very detailed study is conducted. 


The approximations, in addition to the standard ones (ref 10) are: 
1) Forces, moments, and blade motions are quasi-static linear functions 
of the three aerodynamic variables aw,ae , anday, and the two control 
mpuussSo, and OB,y, . Implicit in this assumption is: 
a) the neglect of unsteady aerodynamics, or aerodynamic lags, such as 
the familiar downwash lag, and 
b) that rotor plane dynamics can be described by additional, redundant, 
equations in the helicopter variables. 
2) There is considered to be a plane of longitudinal symmetry, and the 
motions in this plane are uncoupled from the lateral-directioral motions. 
3) The gusts are modeled as moving air mass fronts along ody axes. 
4) No attempt is made to simulate gust alleviation, gradual gust 
penetration, or gust spatial distribution. Gust alleviation, as it is 


used in this work, is defined in Appendix 1V. 





The equations are derived according to the procedure in Refererce (11). 
The aerodynamic forces in the form of Taylor series, expanded in the three 
aerodynamic variables along the X, Z, and about the Y axis are summed with 
the control and inertial forces and moments in those directions and divided 
by the mass and the moment of inertia respectively. The terms are 
linearized for small perturbations and the sum set equal to the trim, or 


initial values. This leaves equations of the form: 
- AU - Xudu-Xwaw +O. AS + QbO = Xz, dB, + X90, 80, 
ae) — Se — WN - By Awe = VU.O = Cremer, + aap d8, 
-Mysa -Mwow —-M&S0 = Ma, dB, + Ma, dd. 


The stability derivatives are dimensional derivatives with the units of 
1/sec for the force derivatives and 1/sec-ft for the moment derivatives. 


They are defined as follows: 


_ 1M 
Ky- & oh Zr k MO * TO 


Oc) aC) 


where the partial derivatives indicate that four of the five variables 
are held constant while the derivative is taken respect to the fifth. 
This does not imply, however, that any dependent variable such as blade 
flapping, as it is treated in this work, is held constant during the 


process. 





The equations describing tne motion of the rotor plane (blade 
flapping) are, consistant with the assumptions made, independent 


redundant equations of the form: 


O,= Wow + Xu AO x OO) ay a Shy OB, 4, AQ. 
Ow Coe OU By 08, 


= Ale AW + Ode - Kd ae Odo -+- he SB a Ji (SO 
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This work pertains to the gust response of aircraft and so a 


Suitable way of describing the gusts, especially for use on the analog 
eomputer and for ease of physical understanding, is used. The aerodjnamic 


weriavle, AW, iS uged tO represent aircraft vertical velocity with respect 


to the local air mass such that 


angle of attack 
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A gust, WE is modeled as the velocity of the local air mass with respect 
memuicstaclOnary atmospnere, and the aircraft velocity with respect to 
mie stationary atmosphere is “We. The sum is the aircraft velocity with 


respect to the local air mass such that w = (avg + awa). 

@ 
This substitution is made into equation (1), noting that the w term in 
me 7 TOrce equation is an inertial acceleration of a body, not an 


acceleration of the air mass, and should really be Tq . VALT of themes 


shown in the figure which follows. 


4 =7 =-WWg 


AK 


qe SS 


MW 
ATF \F\ 
wi \ \\\ 


Wa Body ARIS ve 3 7 pf y . 


The equations can now be written in a form convenient fer pultifniee on 


the analog computer because the gusts are considered as separate inputs, 
or as forcing functions. They can be apremgellas foll@is Gila paverred i214. 


the computer in a standard mammer. 





2 DW = : ‘se i es a 
2 AWo — > « = Lowy AwWy 
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The left hand side represents the accelerations on the aircraft due to 
the right hand side. The analog computer diagram is presented as 
figure (1). The variables printed on the diagram are really "computer 
variables" because they are normalized to maximum values. The analog 
computer responses are labelled in dimensional units and all, except 
where noted, were made with a 73 ft/sec gust input. The values for 
all of the potentiometers were derived from the data from which the 


derivatives of Appendix I were taken. 


It should be noted here that an analysis of this form differs 
from that where the transient response of the aircraft due to a control 
pulse input is studied. A control pulse input will cause different 
relative accelerations, along and about the three axes, from a gust 
disturbance, resulting in different aircraft motions. This can be 
seen as the difference between the numerators of the transfer functicns 
fOr motions resulting from a collective input as opposed to a vertical 


esc Input. 


Ad . Me, (s+T). d6 —. sw 
AB. A(s) D(sq) A(s) 


The character of the motions are the same, but the actual responses are 
dependent on the numerators. When dealing with flying qualities the 
actual response is often of more importance than the character of the 


moulon. 


Much of this analysis will consist of the interpretation of 
Piston cCompucer runs, in Contrast to the solution of transfer functions 
or mode ratios and mode shapes. The analog computer approach is taken 
in order to get a grasp on the basic nature of the responses and because 
of the ease in the analysis. As such it will be physically intuitive to 


keep in mind the simplified computer diagram which follows. 
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In this way one can see directly all of the parameters maxing up the 
airplane accelerations and, of importance in interpreting the modes 
of motion, can see the number of integrations in each loop. It is 
eonvenient to think in terms of electrical signals, adding and sub- 
tracting from one another and being integrated. In the second part 
of the thesis various feedbacks are investigated. Each is merely 

a combination of the helicopter variables summed together and called 
Beeeteedback Signal. This is shown in the full analogs computer di- 
agram of figure (1)where all of the different feedbacks are shown 
@uetvne right side of the figure. Selection of one or more of the 
feedbacks is made by throwing the proper switcn, which sums together 
a particular combination of all the variables and amplifies, attenuates, 
Gr integrates it according to the method chosen. The Signal is 

then put into either the Bj, Alp oF Oo A|p amplifiers on tiewene 


Sieune diapram as a control input. 


Approximations to the equations of motion will be used frequently 
mer Case Of analysis. It will be shown in the first part that the 
approximations used are valid for the purposes of this work. Tnis will 
be shown by comparing results of successively brcader approximaticns, 
and, at the same time, will bring out the major differences between 


helicopter and airplane gust response characteristics. 


It is necessary to have a standard of comparison fcr a study such 
as this. Unfortunately, there is no known single parameter that en- 
compasses all of the important flying qualities parameters. Of even 
itere Concern 1s the lack of knowledge as to what the important fiyane 
Qualities are. For the purposes here, however, rotational and normal 
acceleration magnitudes, as well as dynamic stability characteristics, 
are used aS a general index of comparison. The acceleration magnitudes 
wee WwpOortant Lor structural and ride quality coOnsidemetions andere 


dynamic stability is important for control and flying qualities. 





PART I 


ANALYTICAL AND QUALITATIVE COMPARISON OF AIRPLANE AND HELICOPTSe SUST 
RESPONSE CHARACTERISTICS. 


In this part different approximations to the full three degrees of 
freedom of motion of an aircraft are considered. The value of each 
approximation is shown and the important points it reveals are discussed. 
As each approximation introduces stability derivatives into the equations 
tjmey will be simply derived for both the helicopter and the airplane. 
Typical values for the derivatives are taken from data presented in the 


e@enaices, and in the reterences. 


in erder to gain an intuitive understanding of the various factors 
affecting an aircraft's response to a vertical gust, it is helpful to 
Bearo as simply as possible. Consider an aircraft, or any body with scme 
Jift curve slope represented by the stability derivative Zw, restrained 
from moving in any direction - bolted down in a wind tunnel for example. 
If the body is exposed to a vertical component of air velocity, ce then 
there is an aerodynamic foree created on the body in the Z direction, and 
it .Could be measured by the wind tunnel balance. The force is 

Z = SF BW4 


Divided by the weight gives the normal load factor, 
Ae me 4 
oa S 


If high frequency and, in the sense used here, higher order dynamics 
are neglected (e.g. blade dynamics, structural response), then this is the 
initial vertical gust response of an acrodynamic body. The higher order 
dynamic response is referred to as gust alleviation and cifectivel) er: 
SO as to reduce the magnitude of the initial=response. Trieetoud tae: 


transfer function may be thought of as 


LO 





(a) 01 = — (Geo) 


where G(s) is a dynamic transfer function which, for very low frequency 
inputs can be approximated by unity. (Appendix IV has a further discussion 
of gust alleviation.) This approximation then becomes a static approximation. 
It seems obvious that the derivative, Zw, will be important in the gust 


mesponse Characteristics of aircraft. 


For a conventional airplane it can be shown (reference 11) that 


w4 —_-_ — A Ow = Wie Ow 
iad ene La C 


airplane 


where w= /psz is airplane relative density and a,, is the slope 

qeeune Litt curve for the whole airplane. Zw is proportional to airspeed 
baeeinversely proportional to size and relative density. Small, light, 
high aspect ratio airplanes flying at low altitude and high speed have 
large Zw's (order of-2) and heavy, high wing loading fighters flying 
slowly in approach have low Zw's (order of-0.2). 


Hees voeeaciicopeer, Zi 1s markedly different. If all of tue laa 


Seneirdered to be from the thrust of the rotor, 
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In terms of the rotor parameters Zw is proportional to 
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tic indvecd velocity parameters 


On go) oe of 
Oleoane OWeng OCT 


are taken from reference (12) wnere they are shown to approach UNG” 
Zero, and zero, respectively, as forward speed increases and induced 

effects become negligible. Using these high speed approximations and 
mMoting that, in level flight the second term will be negligible with 

mespecy tO the first, then 


Cc = 
ocr — —_ for high speed. 


dw 4 C22 


Dimensionalizing yields 


acy 9%) 
(RR) (SY) 


helicopter (high speed) 
which is similar in form to that of the airplane, but is independent of 
airspeed. For a helicopter Zw approaches a constant value, for a given 
blade loading and RPM, at some high enough forward speed that induced 
effects are not important. It can be seen in the data of Appendix I 
that the speed range considered herein is not yet high enough to warrant 
this approximation because Zw for each helicopter is still increasire 
even at 180 kts. It is also interesting to note that the helicopter 
fepresented by the circle and the square are the same helicopter where 
the larger Zw corresponds to an empty weight condition and the lower Zw 
to a heavily loaded condition. The low "triangle" values are for another 


similar helicopter which is at its maximum gross weight. 





In order to contrast airplane and helicopter Zw's it must be noted 
that Zw is highly dependent on the flight condition for a given airplane, 
whereas for the helicopter Zw approaches a constant with speed. Jhe 
@esien ilight condition thus enters the picture when a comparison is made. 
It is clear that at some speed and altitude (dynamic pressure) the heli- 
copter and airplane will have the same Zw. Figure (2) indicates that, 
for six dissimilar aircraft for which data was readily available (ref. 10), 
Miesspeed waecre this happens may be near the design condition for the 
Particular airplane. The values of Zw for all six aircraft near their 
design flight condition were within the range of data for four similar 
helicopters at high forward speed. For purposes of comparison, it is thus 
considered that an airplane, designed to fly in the speed range of interest 
iemclis work, is likely to have a value of Zw similar to a helicopter flying 
av vhe same high speed. Typical values are on the order of unity. It is 
fame concluded that the static load factor for an airplane and a helicopter 


at the same high speed is similar. 
The next approximation is a single degree of freedom approximation in 
fren Only vertical motion is allowed, the other two degrees (pitching and 


horizontal translation) being restrained. The body can be thousht of as 


ieaislatving On a vertical rod thru its center of gravity. In this case 
LU NOVO) 


so, of the equations of motion, the only one that applies is 


Wy : eo OQWy = Zw ae 





The aircraft response to a steady vertical gust is 


“bw 
Wa = “Wall © "floating" 


ene 
No = - £w Awe Cc normal acceleration 
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This simple approximation shows that the time constant of normal 
acceleration is also dependent on Zw. In the figure are shown typical 


responses for different Zw's. 






DE Larqe Teil 


Small ~w 


+ wae 


The large Zw causes a higher peak load factor, but a faster tine 
constant in reducing the load to zero. Statically, then, a small Bw 


is desirable, but dynamically the larger Zw is seem tom¥e betmer bee 


cause of the increased rapidity with which the gust ifduced acceleration 


are reduced to zero. 


Ly 





The approximation obviates the fact that a lifting surface, airplane, 
or helicopter will "float" up in response to a vertical gust until the 
aircraft's vertical velocity is equal to that of the gust and the engle of 
epcock 1S thus zero. The longér it takes for this to happen, the longer 
wee une aircrait experience a normal acceleration, as well @s ovher 


i] 


forces and moments due to a non-zero angle of attack, (e.g. M,ox, X,o%). 


From the first two simple approximations it is seen that Zw, tire 
mee, curve Slope of the aircraft, controls the load factor, and that 
euwenavic direct Jift or thrust control could be used to either decrease 
the static load or decrease the time constant of the response, but could 


not accomplish both objectives at the same time. 

The next step is to consider two degrees of freedom, pitching and 
heaving, at constant forward speed. This is the so called short period 
approximation of the conventional airplane, but really is an approximation 
to the motions ensuing before airspeed has a chance to change appreciabdly). 
whe approximation is exact for the time it takes for three integrations, 
if Xw = 0. (This can be seen on the simplified analog diagram in the 
Bescription of Analysis Section by tracing a signal entering as Ng and 
iiadly causing a AUg ). 

The equations that govern these motions are 

a We — 4s Aude Ue SQ = ow 5 


sO ~ M@6AG - Mw dbWy = Mu Awa 


ema wne characteristic equation is 


(s-w)(s-vs) -Me = 6 


“ 





It was found that Zw is likely to be of the same magnitude for the 
helicopter and the airplane, and now Mg and M, must be investigated 


for contrast. The airplane will be considered first. 


For an airplane the moment derivatives depend heavily on the 
horizontal tail. The only contribution to Mg usually accounted for 
in estimates for "conventional" airplanes comes from the change in 
tail lift with tail local angle of attack as the airplane pitches. 


One way of writing the expression for Mé eS: 


—_——S— - 


Me = de Sa BH NS, PU Rce| 
Ty du, 96 A Ty “eal 


Mé = Li UO: Gare 
A Ky MA 


Lt airplane 


where U, =U.(tail) (i.e. the absence of power effects) 


and 


V. Wh ae le °) Ci = Oy NE 


So 


The magnitude of Mg for the airplane thus depends on the tail size and 
the square of the tail length, and is proportional to forward speed. 
Typical values in the speed range of interest are between -1.0 and -2.0, 


as taken from the data in reference (10). 
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An expression for M, is easily written from a force diagram by 


summing all the moments about the center of gravity 
2M=oO= (Xce - Kac) Ls ei M ¢ a Genie 


and taking the derivative 


OM = (Xce 7 Xac) 4° Ow + q=c Con ck g ~ f[decv) de so 
Ok 
In the absence of power effects q = q (tail) and 


ey 4 
Maa = Gi (Kee = Kec) Oy + Cun oe 5 -_ Ce 


2Ky A 


For the "conventional" airplane, Me is thus proportional to dynamic 
mresosure and depends on center of pravity position, tail lengtn, tail 
Size and fuselage characteristics. The tail term is usually the dominant 
one in the expression such that My is large and negative to provide 


ample static stability. 


The "conventional" helicopter is not at all like the airplane in 
regards toM,and Mg. For the conventional single rotor helicopter 
With no hinge offset or horizontal tail, and neglecting the fuselage, 
Mg arises from the lag in the tilt of the resultant thrust vector as 
ine sSshatt pitches. For these helicopters the lag of the thrust 
vector is approximately equal to the lag of the tip path plane of the 
rotor (reference 1). Mg arises from the increased drag on the rotor 
blades (H force) as angle of attack is increased. Both of the 
dérivatives depend on hub height, and are independent of center of 
gravity position. For the hub above the center of gravity io is 


negative (stable damping term), and My, is positive (static instability). 


if 





The helicopters capable of flying in the 150 kts range, and especially 
the large, heavy ones of this study, are very unlike the "“convertional" 
helicopter. These newer helicopters have main rotor flapping hinge offsets, 
horizontal tails, and large, bulky fuselages. The reason for the hinge 
offset is for greater control power, and the tail becomes necessary (and 
effective) as speed increases. This work is concerned with only these 
latter types of helicopters and all reference will be to the helicopters 


mapresented by the data in Appendix [fI. 


From the force balance diagram in Appendix (II) comes 


SM=0O = TLR -h(Bs +i.) | + Hf rh (Bseis) | 
i (1) 


ay cht? a, x My i he oi 


; : : », 
When the small angle approximation, H = Ta is made, the derivative 


Taken with respect to 0, and it is noted that 





then 
OM. AT da 4 cbMsY dar — Me 4 Sk a, 
© eC a © Le). 


From the relation a’ =Cu/ey , 2b Can be “shown thav (reference 1) 





Cee Ca Se se -~ A dal 
QO 


31 GU a 
a6 Swe Cry, OU 
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Qo Cr 


As fe increases, as with heavier faster helicopters with the 


on/ 
rotor plane tilted farther forward, it is seen that OO decreases. 


It has been found characteristic of the helicopters investigated that 
pitemnern results in a small contribution to the pitch damping at 150 kite, 


' 
The trend of O8f- with speed for the data nelicopters is shown in 
de 


figure (3). For the 24000 1b helicopter, for instance, this term (which 
makes up all of the pitch damping in the "conventional" helicopter) 
meaerrouves less than five percent to the total damping. 

From this information it is seen to be reasonable to neglect the tilt 

of the rotor resultant force vector in estimating the pitch damping 


at high speeds. 


The time constant of the rotor tip path plane tilting response to a 


aero Lint rate-is 


204 age a - + Oa, 
Oe sh. (1- = ) OU 


which, for the speed range and helicopters of interest can be approximated 


toy eg: Including these assumptions results in 


_ Seb. Q’ _ (V2) Whe ease: Ors y 


Mo 
¥Ty Ly 


tas 
where y= a / © the ratio of the dynamic pressures in and out cf the 


rovor slipstream at the tail. 
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For the data helicopters the table below shows the approximate 


breakdown for the contribution to Mg at 150 ctaes:. 





Oo 
HELICOPTER (C+/0) Me TAIL TERM HINGE OFFSET TERM 
2000 fwd gyi BiG 15% 20% 
35000 aft 3.0 al 60% 35% 
42000 fwd emis . 38 Loz. Lo% 
18000 aft DOT, 53 50% 35% 


Typical values for the pitch damping are seen to be between -0.5 and -1,0 
at 150 kts. It is clear that the need for a tail increases with Bo /er 


ave - | = 
(decreasing OC ), and the data is in agreement with tre prediction of 
reference (1) which says that the tilt of the thrust vector contribution to 


tie pitch damping should be zero for 


Of) = oe 


Peoetal. Size increases the pitch damping approaches that of the airplame 


with the same tail size. 


An expression for My results from taking the derivative of equation I (1) 


with respect to & . 


OM = iy Oo ee -4- L (Bis ¥is) -- a'{ TBs ri] 
OW 4 OK 


: 4 + [Be ris) he elo OF thy } oA a ot 
2 O 





Making the approximation that, for any helicopter => An (Sis ris) . 


OCT, 
using the high speed, approximate expression for Oy) derived for Zw 


earlier, and making the approximation that, in the absence of induced effects 


Oa aes 2 aPN — : 
ame \ 50 te 


eoen, the expression for M,, becomes 


Yo aon ae rao: » chveeier 
X 





se Ky GOR) CY) Ty 
3 
+ a Y Cw ss Cur ec « 
3 
where 





TR? 


éE = g' a - se (eet. | 


AlY of the terms except Cui and, pessibly. the secondeciesane 
destabilizing. As in the airplane, the tail should be large enough to 
evercome the instability of the rest of the aircraft. Unlike the airpiage. 
lgevever, there are other restrictions, such as control to trim al difiereue 
power settings, that govern the size of the tail and have kept it relatively 
small on all present helicopters (relative to the airplane of the same gross 
weight). It can be seen from the expression above that there is little 


chance for M pg to be negative without an extremel;, lange tail. 


el 





It should also be rememoered that the fuselage is characteristically large 
and bulky for these helicopters, resulting in a much larger fuselage 
contribution to the static instability as forward speed increases than for 
the airplane of corresponding gross weight. The argument, from a staodilit; 
jeune: Of View, is for a much larger tail, but as said, other factors are 
likely to keep the tail size smaller than on the corresponding airplane. 
The tilt of the resultant thrust vector term and the hinge offset term are 
beth destabilizing, and both increase with a power of the forward velocity 
as can be seen from figure (4). The second term in the expression 

@en be stabilizing or destabilizing, depending on the orientation of the 
resultant rotor force vector, R, with respect to the center of gravity 
memos trimmed tilieht condition of interest. The term is stabilizing if 
ie extension of the vector, R, is behind the center of gravity, € is 


negative, as in the figure. 


ANF f2 SHARE ORS. | op ypP 
» 


4, 


te 
V 


i} 


yo fe 
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In order for this to be the case, the net moment in tHe trimmed flignt 
condition contributed by the flapping of the blades about the offset 
hinges, the fuselage, and the tail must be positive. This means that the 
tail should be at negative incidence and thus carry a download to increase 
the size of this stabilizing term. <A dcwnload is also desirable to in- 
crease static speed stability, Mu, but other factors such as dynamic 
Staoility, and control required versus trim speed and power determine the 
incidence of the tail. It should be noted that, for the "conventional" 
helicopter, neglecting the fuselage, this term is identically, zero because 
Waepmorce veecuor, RK, must pass thru the center of gravity in trim. 

For helicopters with hinge offset, tails,and large fuselages, however, this 
tem 18, in general, non-zero and increases with forward speed. It is also 
Seem tha. the center of gravity position plays a smaller role in the static 
stability of the helicopter than for the eirplane. The C.G. position, by, . 
only occurs in the € term and its importance is obscured and diminished 


as all the other terms in equation I (2) get more important. 


The actual value of My is very dependent on the particular configuration 
as each term in the expression can be of the same order of magnitude. It is 
probable, however, that My will be positive or unstable for present heli- 
eopters in the high speed range. The average value for the data relicopters 
Prel5O KES 18s about +2.0. 


Returning now to the characteristic equation for the constant speed, 
AU= 0, approximation, if "typical" poles are plotted on the complex plazie, 
and the roots traced as Mx is varied to "typical" values for the airplane 


and helicopter, the differences in characteristics become apparent. 
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Helicookec v airplane 


For the helicopter, instability occurs when 


(Mx a Z.. Me) pe 


which, for "typical" values of Zir and M@ given earlier, is when Ny 

mee voerween +O.5 and +1.0. This is below the values for all but one of (Cie 
data helicopters, and the effect can be clearly seen in the analog 
computer responses for the different helicopters in figure (5). In all 
but the most stable one the motions are dominated by a rapid divergence 

in pitching - a consequence of angular acceleration induced by 


positive My. 


For the airplane ,the roots are of short period and well damped, 
iying in the region indicated on the figure. This mears that the airplane 
will accelerate to the gust velocity very quickly, and the moments and 
forces due to angle of attack will be quickly dissipated after a 


disturbance by a steady vertical gust. 


oO), 





A further approximation to the equations of motion to adequately 
describe airplane accelerations in response to steady gusts is not needed. 
The airplane typically exhibits a phugoid, or long period motion, in AU 
and A@ while 6W= 0 (in the body axis system). The resulting accelerations 
on the airplane are small compared to the ones in the "snort period” 
response. For the helicopter, however, the pitching and angle of attack 
responses are influenced by the changes in airspeed because of the neture 
of the derivative My, the forward speed static stavility. For the 
"conventional" airplane, Mu = 0, but for the conventional, es well as the 
high speed large helicopter, Mu is in general non-zero. An expression 


for Mu can be written as 





OM re KT Ia + eLMsQ- Oa 
U 





Boe Cv Ov 
+ TR? ES ey _ Ve a 2 ie 
4 Ou v; Fuselage “ ou wy: 


OCur OC 
where ( Su s+) (fuselage) and ( 4 ) (tail) are influenced b: induced 





eizecus or changes of the rotor slipstream with speed. It can be seen 
that the size and magnitude may depend largely on the incidence of the 
Merizonval tail, if the tail is large. Positive Mu is statically stable, 
but, as will be shown later, positive Mu leads to dynamic instability. 

All of the terms in the expression vary scmewhat irregularil, with 
forward speed because of the induced effects of the slipstream, even at 
Speeds as high as 150 kts for the data helicopters. With a large enough 
tail, at high speed, the sign of the term can be assured by adjusting 

Waer vail incidence, but, with a small tail, the value of Wm is vemy 
dependent on the particular configuration and flignt condition (e.g. level 
or non-level flight). Of the data helicopters, the two that exhivited the 
highest and lowest values of Mu both carried tail down loads, while t8estvo 
that exhibited little chanre Witl] speed catmied a sig postive or wero 
load at 150 kts. The smallest #u's were for the smallest Rela caver saa 


Gne smailiest Lai Ll. 


nN 
Wa 





An approximation that includes only this particular effect of the 
airspeed changes is presented as follows: 
It is assumed that Xu, Xw << g and Zu = QO. 


The characteristic determinant can then be written as 


S O 5 
O S- tw -U25 = Go 


“Myo -My = S(s-Mo} 


piemtMe characteristic equation, in root locus forn, 
5° ( AU=O AP PRO K ROdTS ) “a aMu (s- %) = O 
For a positive Mu, which is taken from the data as typical for these kinds 


Seenelicopters, the root loci for increasing Mu will look approximately like 


one of the following 





AW=O APPROX (0) \ Au=0 AfPRIx (a) \ 


STABLE UN STABLE 





From the figure it can be seen that, regardless of the 4U= 0 root 
Hocations, the helicopter will always exhibit two real cconvergences if 
fies positive. These roots can loosely be thought of as the convergences 
of angle of attack and pitching rate, the larger one being the faster angle 
of attack response. The inclusion of the other 4u related derivatives is 
shown in the analog computer responses of figure (6), where 
a) is for a helicopter with the "typical” derivatives (except for Zw and Mw) 
b) is for the lowest values of all the derivatives (except Zw ard My) and 
c) for their highest values (except for Zw and My). 
It is seen that there is very little difference in tne acceleration 
responses, Na and6@, for the first few (+10) seconds throughcut the 
range of all these derivatives. The reason for this is that the rate of 
convergence in angle of attack depends primarily on Zw, and the accelerations 
are primarily dependent on the angle of attack. For the first few seconds 
the acceleration response is almost independent of anything except Mar and Zw 
and, as such, should be accurately approximated by the two degree of freedom 
model in which AX= 0. This corresponds to the fast convergent root which 


represents the A& motions before coupling witn AW motions. 


The result of all of the discussion above is the point that there should 
be little difference in the first few seconds of airplane and helicopter 
normal acceleration responses, regardless of the stability or instability 
of the two degree or three degree of freedom models. The static and the oxe 
degree of freedom models show, with the assumption that Zw is of the same 
Magnitude, that the responses will be identical. Inclusion of pitching as a 
degree of freedom only changes the time constant of the normal acceleration 
convergence, making it longer with positive My and shorter with negative ty. 
If Mx is large and positive the convergence becomes a divergence and 
anclusion of airspeed as a degree of freedom is necessary to describe tke 
motion, but when this is done (for positive Mu) it is seen that a convergence 
reappears on the complex plane. Physical interpretation of this will be 


given in the following qualitative discussion. 





Regardless of root locations or characteristic modes of motion, there 
are some fundamental differences between airplane and helicopter pwust 
responses that do mot show up on the complex plame. The biggest, mest 
@evi1ous, and most Tundamental is”tne fact that with different signs of 
epee ecopLerseema alr planeemprich in dilferent directions in respomse 
to a vertical gust. If My for an aircraft were zero, its initial response 
to a vertical gust would be pure heave - no pitching at all - and a very 
slow and small response in Au whose sign depends on Xw. (Xw can be 
positive or negative but is always small.) For the airplane with negative 
Peel eeeresponse 15 2 nose down pitching along with the vertical acceleravicn, 
This causes an increase in airspeed as the important 4U response (for 
Xwe< eg) due to the eae) term in the X force equation. The airplane then 
SibecremmnLO MesmelWaracterisuic long period, or phugoid, oscillation. Gee 
opposite initial response occurs for the helicopter with positive hy. 

The helicopter pitches up, and the important 6uy response is a lcss of 
cumespeed. in addition, if Mu is positive, as it is with all of the data 
helicopters, the AU response decreases the pitching acceleration and tends 
to cause an oscillatory mode if the effect is strong enough. The following 
diagram illustrates the differences between the two aircraft by showing what 
an Observer, flying straight and level at a constant airspeed, would see 


when the two aircraft encountered gusts. 
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When the helicopter oscillatory mode is of a long enough period, the 
Mejor part Of the acceleration responses rill be convergent, like the air- 
feene, and will take place in the first fewsseconds. A pilotew#euld 
lecovevly now notice the differences in the normal acceleration respoemse 
for the two aircraft. In figure (7) are snown analog computer responses 
fore une same aircrait but with 
a) positive My, 

b) Me= 0, and 
c) negative My. 
Both the constant speed approximation and the three degrees of freedom 


Hesponses are shown, and they illustrate the observations made above. 


Wea conclisyon to this part of the work, it is possible to summarize 
the findings qualitatively as the differences a pilot would eucounter in 
tie two types of aircrarIt as a result of a gust disturbance. First, he would 
imeeopabky MOL nOLice the difference in normal acceleration, either the 
maximum load or the first few seconds of the dynamic response. Second, and 
M@sveimporvant, he would notice the different direction of rotational 
eee rariOn, especially if he is far ahead of the center of gravity. In 
the helicopter he will experience a local "g" response greater than that of 
the airplane because the helicopter pitches nose up. Third, he will wotice 
a smaller angular damping and a tendency for the helicopter to continue to 
rotate until the static speed stability causes restoring moments, usually 
Resulting in a long period unstable oscillation. Finally, because the pilot 
meeoole to respond with his controls, he will notice the difference in 
eovurol response required of him in order to dissipate the accelerations. 
In the airplane the pilot pulls his stick back to keep his nose up in 
mesponse CO a vertical gust, and in the helicopter he must push fom@ard 


on the cyclic control to keep his nose down. 
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PART IT 


ROTOR FLAPPING FEEDBACK CONTROL OF GUST INDUCED HELICOPTER DYNAMICS 


This part of the thesis pertains to the elementary stabdilization of 
me heavy, high speed, single rotor helicopter which was characterized 
in the first part of the thesis. Tne stabilization schemes that are in- 
vestigated include a conventional rate and attitude system which is 
compared and contrasted with various types of rotor flapping, tilting 
moment, and thrust feedbacks. The impetus is to determine whether or 
not the latter are adequate to provide desirable gust response qualities 


mom a helicopter. 


It is often considered desirable to reduce the amount of rotor 
flapping that occurs due to a gust disturbance. This is essentially 
an attempt at devising an “inherently stable" rotor with which to 
stabilize the helicopter. The reduced flapping is generally favorable 
to blade motion stability, but the extent to which this “imoroves" the 
@ymanies of the helicopter is not given proper attentior. Delta Three 
and Oehmichen pitch flap coupling, and variations of these, tne Lockheed 
gyro system, and a hub moment sensing system are examples Of eacse 
techniques. In most of the literature, however, only the “fixed shaft" 
dynamics of the rotor, equipped with the particular device, is studied. 
The consequence of this is that the actual effect on a helicopter's 
Mena Ggualities 14S not known. In this part of the thesis assim yea 
investigation of the effects on the helicopter as would Se encountered 


in actual flight are investigated. 





Influence of Rotor Flapping on Stability Derivatives 


The fundamental problem in trying to determire how rotor flapping 
influences helicopter stability is to determine how much of the totel 
forces and moments on the nelicopter is attributable to flapping. An 
assumption often made with the low speed conventional helicopter is 
treo the rotor resultant force vector, R, is perpendicular to the 
rotor tip path plane, and that the motions of the two coincide after a 
Gisturvance. This assumption leads to the erroneous conclusion that 
the position, or rate of change of position, of the rotor tip path 
flane due to a disturbance or control input is more important in the 
Meese and Control cf the helicopter than it reall: is. Stability 
only is considered in this work, but the same arguments as will be 


presented can be used for synthesis of helicopter control s:stems. 


Precisely speaking, helicopter accelerations only partial, depend 
on rotor flapping, and as power, speed, and size increase, these other 
factors become of prime importance. Because of the way in whicn rotor 
flapping is treated in this work, i. e. as a dependent variable whose 
magnitude depends only on the helicopter variables, the influence cf 
flapping on the stability derivatives is implicit. In order to expose 
the flapping dependent contributions to the derivatives, it is necessars 
to consider the flapping angles as "quasi-independent" variables and to 
take partial derivatives holding them constant. A moment derivative will 


thus Jook like 


OM “= OM. ae om On + OM Obi + om othe 
30) OG) Wf dC Ob OQ) dds} 90 
0,0. b, hs bo, AG, Gib, 





Expanding, because a' = a' (a,, a,, b, ) results in 


o 3 


ees 4 24.) My] dai 4 dM 
OC) OC) C\{ da’ aoe 04, F 


Ab, as b, Qo 


rom) oC) da'| Ob: db, Pla 


p OOo, OM. a + OM 


OC){ da 20. DOs] - 


Bach term in the equation has a special significance, whtch is listed below: 


OM Change of moment with ( ) independent of 
OC) rotor flapping. Fuselage, tail, tilt of the 
a, b, Qo 


R vector independentl;, of the TPP, andy vugrae: 


changes make up this term. 


\ 

OM OO, That component of the tilt of Wie R vecvor 
\ 

OU Od, (H force) dependent on the longitudiral 


Cilgeaor-the Tr. 


\ { ¥ 7 Vix 
OM O6, °) OM Ou! That part of the tilt of the R vector (E force) 
; , 
Ou Odo OU Clo dependent on coning and lateral flapping angeles. 
OM Divect hub moment due to hinge offset (or 
eee \ 
OQ { flapping restraint, or any oter darect Guage 


of moment due to the wechanical aspect of 
flapping as opposed to the aerod;'Nawic forces 


caused by flapnpi y ) 


OM 5 om Lore . (: ll nelWtiecanmser ver 
Ob, Othe cone tt t | reertir) 
, 





The term, oa" | can be thought of as a measure of how closely the 
Oa, 

mower resullanly force, R, and the perpendicular to the tip path piawve, 
fMeeeecoierde. it is like a correlation function in that when tke 
derivative is unity the R vector is perpendicular to the TPP ard when 

the derivative is zero the orientation of the R vector is independent 
faerie IPP, The terms Oa, and oa/ ean be thought of as dihedral 

Qo Ob, 


effect due to the coning and longitudinal blade flapping velociiioaye 


The following figure illustrates this. 
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Tt is seen that the longitudinal flapping velocity (lateral pula ee, Dw 
meauceswunesdinearal Coning angie etiect, thus decreasing the bdaGigard 
mie Of the R vector. It will be shown analytically in the coming 


discussion that the two effects almost cancel. 


The important distinction to be made in this analysis is the 
se e ] e e 
partial independence of a and a, , and the important result is 


"by how much?". The angle a’ can be written as 
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From the expressions above can be found the derivatives 


OO. *) 0! 9 Qa 
OO, On, Ob: 


They are, after much algebra: 
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Hoarcuhermore, from the expressions derived in Part I for the moment 
balance come the following: 


OM] _ wT OM) 2. ebM, 2’ 
Oa’ , Oth, a! 2 
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Substitution of all of the expressions derived above into ecuation II (1) 
yields 


OM. = aM 4 dts yal + QO & + 0,-0.4] 
dC) ad) OU) Gey \ * 
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The last term in equation II (3) gives the net effect of conirg ard 


lateral flapping, and can be expanded as 





ATar| qQ. dO. \ ISSR Pe + 
12QCr AC) (y+ hp") 9 


ayo. | ieee Vela 
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When the variable to be inserted in the (_ ) isQ, the expression is 


approximately 


Gee 
fae? Cy 


but for any other variable inserted in the {( ), the expression 


reduces to 
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item Lovcmabe stial! in any condition. . The effect of lateral and 
eollective flapping on the moment derivatives is thus seen to be small, 
e-pecwo) lyawben 1t is memembered fyom Part I that, for the high speed 
helicopter with a relatively large tail, mos& of the pitch damping 
comes from the tail. The last term in equation II (3) is considered 
zero and means that coning and lateral flapping cancel each other's 


influence on the longitudinal dynamics of the helicopter. 


i wae Gete Liao bine moment derivatives can be written as the 
MimotetyO vermis, Chose depending on longitudinal flapping and tnose 
not depending on any flapping. The flapping dependent contribution 
ean further be classified as that part due to flapping hinge offset 
Wiemeiueearo due to Lhe coincidence of R vector tilt with TPP tilt. 
Peeelabter as the term whose magnitude is often overestimated by the 
RLTPP assumption, and whose importance will be investigated in the 


eoming discussion. 


It has been said that the coincidence of the R vector and the 


SPP movlons 1S measured by 


doo = 4 yg OE PR 4 ma - Ou it (2¢) 
O4\ Ge, | 2 
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At very low inflows and forward velocities it is easy to see that the 
RLTPP approximation is fairly good. As speed, power, and inflow in- 
Seeace, NOowever,, as they do with the newer high speed helicopters, re 
Sign and magnitude of the sum within the bracket is not obvious. It 
bem be said, however, that in a high power climb condition » and Qo 


fae seun be Llarse enoush to assure that oa, me 4 
\ 


whereas in a low power descent, autorotation for instance, A and Oe 


wilj. be small and do 
Mya, 7 4 


In level flight it is necessary to expand tne terms and make 
approximations to the expression in order to gain insight as to 
what are the important trends and contributions in terms of more 
Hundamental parameters. The obvious sensitivity of the expression 
wemertin Conditions should instill caution when considering approximations. 
Usually good approximations are found to alter the magnitude ard even the 
Sign of the whole expression because the terms are all about the same size. 


Substituting into the bracketed part of equation II (2c) for a Q, gives 


Da 


L 2 tye - 064° | a eP 
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in which the approximation 


[i- 5p | |) 


is made (normally <5%, error), and into which is substituted 


as - es \ -> | 
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Finally, the approximation 


eee | 


Grace 


is made (normally <5% error) and the resulting expression substituted 


back into equation II (2c) to give 





7 Nee ee 
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The inflow ratio, A iets 


> = Jasen Oper z ik 
Oy (= 


and is obviously negative in level flight and grows with forward speed, 
but the actual magnitude is sensitive to angle of attack and induced 


effects. Because of this the trend for a decrease in Ow « 
) 


with inflow is seen, but magnitudes in level flignt are still not 
obvious. The values exhibited by the data helicopters are presented 
in figure (8) and will serve to characterize the "typical" trends for 
heavy, high speed helicopters. In table II (1) are the values found 
using the exact expression given in eduation TI S ab LDOo ts. 

The sensitivity is seen to be great, for the approximations are medest. 


The approximation does, however, exhibit the right sign «nd trends. 


Piste I (1) 


HELICOPTER x da! Aa! 
7 (Cr/) aa, a, Cores) 
42000 FWD Si Se BESS. 
24000 FWD -1.16 35 4h) 
35000 AFT -1,.26 250 oT 


18000 AFT =? et = 04 aah 
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Theweeyecs Of iis analysis is to find wnat direct effects rotor 
flapping has on the stability derivatives of the helicopter, not to 
prove or disprove the validity of the RiII[PP eapproximatior. Because 
timers iswotften a point of misunderstanding, mowever, it is considered 
worth while to dwell on the errors tnat would be incurred by the use 
of this approximation, These errors are of special concern to those 
who would devise some kind of helicopter stabilization system based 
on a reduction of rotor flapping. It will be shown that serious 


over-estimates of “improvement” are the result of this thinkirg. 


When the RATPP approximation is made, eauation II (1) reduces to 


aM dar + ebtQ"! dar 4 hel dT 4. at 
Ov) aA OC) o( | CO) (fuselagettail) 


when it really should read 


OM . &T dot 4 ebi4s2 da 4 he OF 4 val 
OC ) A) a Ol) OC) O() (fuselage+ tail) 


/ e e e | ® 
where € is computed with a, instead ofa . (The error incurred here 


is indicated by the plots of a' and a, in figure 9). 


oO 





When the term AT da) 1S Of Significance to tie whole expression, 
OC) 

then the approximation can cause errors. It was seen in Part I that the 

Influence of the R vector tilt in Mg is small with high inflows, or 

power settings. The approximation thus causes little change in the 

value of M6 in these conditions. As the power settings are reduced, 

@eecourse, the R vector tilt 1s more important, but the approximation 

becomes better for the same reasons. The derivatives ™ and Piu are the 

ones most affected. In these derivatives the "nT" term is of the sane 

magnitude as the other terms, and the differences in dO. and da: , become 
0) 00) 

large. Figures (3, /+ and 10) show the divergence between the two 

iemeameters, and indicate the magnitude of the errers. Thewstatieestability 


ms over-estimated in both cases, the error increasing with speed. 


iewapproximeyion deesenot appreciaoly effect any of the imporvagt 
force stability derivatives. Flapping has an influence in the X force 
Gerivatives, Xw and Xu, because of the influence of flapping on tke inplare 
H force. These derivatives are of relative unimportance to this work and 
so the effects of flapping are considered to be only on the mcemert 


derivatives. 


Having now seen how flapping influences the forces and moments, and 
thus the stability derivatives, it remains to be snown by how inuch the 
derivatives may be changed by controlling tne amount of flapping - as with 


a tlapping feedback system. 





Longitudinal Flapping 


Consider first an hypothetical perfect feedback system designed to 
keep longitudinal flapping zero in some unspecified manner. If this were 
an infinitely fast feedback loop such that a, 0, then in the expressions 


for the moment derivatives 


Omeenay = Of) _ OA — 


-_ = = —_ Se me, Oy 


OK Ow ‘Oe OO, 


The derivatives are changed by the now zero contributions from flapping 
hinge offset and the R vector tilt due to flapping. The moment derivatives 
are all reduced and the changes produce a different set of helicopter 
Elearacteristuics. it 1s to be noted that this is not a feedback to the 
eyeclie control (which will be taken up later) but rather a computation 
of stability derivatives ignoring any contrioutions from longitudinal 
fijpping. As such, it is now to be shown that the effective result, in 
the high inflow/speed Condmuron, 15 helicopter Stability and Contra 
characteristics similar to the same helicopter without flapping hinge 


offset. When equation II (3) is rewritten in a simple form, 


Maz) BE dw 4 eb} de 4 1 OM 


ie 5 ) 
Ty da, 1 Ov) Iy a) : 
eam) 
pear hay @ offse4 
R vector 
tilt 
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it can be seen that, if the first term in the bracket is small with 


respect to the second, the net effect of 8a, Ay 5 = 0 


is just the loss of the hinge offset term. Indeed, for the data heli- 
copters in the high speed flight condition, where ae is gmair, 
the following table shows the percentage of the first term to the 
total flapping contribution. 


ioe It (2) 
42000 35000 2h000 18000 
FWD AFT FwD AFT 
% (R tilt due to flapving) 


otal flap. cont. ) 35% 208 ie Ot 


The hinge offset is included in the design of the heliccpter primarily; 
to increase the pilot's control power, or the amount of angular acceleration 
he can impart on the helicopter for a given stick deflection. If this 
Meritect servo-system were installed, however, the pilot would lose the 
additional control power provided by the hinge offset, as well as thew 
mertion Of the R vector tilt due to flapping. The pilot would feel 
almost as if he were flying the same helicopter but without hinge offset, 
especially for the last two entries in the table. This particular feed- 
back thus seems to contradict the reasoning for designing the hinge offset 


mothe first place. 





This same type of contradiction will unfortunately arise in any 
feedback system design in which there is a simple mixing of the pilot 
anputs with feedback inputs. For example, consider an aircraft auto- 
pilot which can sense @Q, Q, or re and feed back the signal to sum with 
meLOol Inputs. If Gis fed back, the pilot is limited in how much torcue 
he can apply about the aircraft center of gravity, but any gust induced 
torques are also limited. Similarly, for ) and © feedback, both the 
pilot and the gust disturbance are limited in how fast and how far, 


respectively, they can make the aircraft rotate. 


It has been necessary to limit the scope of this worn, and thus the 
reduced control problem associated with increased stability to gusts is 
not considered. Certainly, compromises will have to be made in any 
control system such that the pilot has the final authority in the systen, 
if he wants it. In this work the pilot has pushed a "hold" button and 
@esires to retain an equilibrium trimmed flight condition, and never 


enters a control input. 


Along the same lines, it is very interesting to note that a Negative 
a —-=B, feedback system (opposite sign of the one just considered), a SO 
called "biased cyclic control system", has been investigated in reference (4) 
in connection with improving the control characteristics of a helicopter. 
With this negative (really positive in normal servo-convention) feedback, 
and with a suitable linkage, it was possible to reduce the control 
sensitivity (maximum pitching rate per unit stick deflection) by increasing 
the pitch damping at constant control power. It is mentioned in that 
reference that the feedback system seems to the pilot to be like adding 
hinge offset, confirming the assertion made in this work. Gust imduced 
loads were not considered in that reference. It is felt that if they had 
been, the system would have been considered undesimasle due to 2nereesge. 
dnostatic instability. There is little feasibility that (Ast parereula 
feedback system, in the high speed flight condition for Nelicoptexs with 
relatively large tails, will improve caijtrol characteris 


1s made up primarily by the tail comfrilmlion, as was shen ig Set 1. 
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At 150 Knots the following changes in the moment derivatives of the 


data helicopters are found due to the perfect, infinitely fast feedback 


loop. 
TABLE II (3) 
42000 FWD 35000 AFT 24000 FWD 18000 AFT 
0,*0 Q,=O Q*#0 a =O a0 G,=O Q,+0 Q,=0 
Mu -—s-. 003 0017 .0018  .0009 002 . 0009 . 0009 . 0003 
Mg -.38 eee ~ 54 . 386 -.767 -.583 8 eee 
Wy. 505 -. O41 7.54 1.98 2aNe -.310 ce son 


The changes in the derivatives change the helicopter characteristic 
equation, and the resulting root locations are shown in figure (11) for 
all of the data helicopters at 150 kts. As can be seen from the changes 
mierooL Locations, improvement in stability is not impressive, especially 


in the aft center of gravity cases. 


The reason the improvements are not impressive is because the fuselage, 
tail, and non-flapping-dependent R vector tilt ( 20/4, small) are large. 
This should be considered as typical of large, heavy, high speed helicopters, 
and will become more so as helicopter configuration approaches that of a 
@ompound helicopter. The truer this is, consequently, the less eitec tie 
will be any type of helicopter stabilization system involving rotor flapping 


feedback control. 


Ny 





For the one case where this feedoack has succeeded in changing 
the sign and magnitude of My enough to put the roots in an oscillatory 
region, it should be noted that the frequency of the mode is on the 
order of an airplane phugoid mode, and does not resemble the airplane 
classical "short period" mode in regard to freauency or damping ratio. 
m@eschange in sign of the effective HM, for both forward C.G. cases 
does change the direction of the initial pitching, however, and would 


result in a decreased "local g" response that a pilot would feel. 


Micmwcemeivcus Presented above 1s for orly one flight condition, 
ao 150 kts, and shows that helicopter stability can be improved, Nowever 
medestly, by reducing longitudinal flapping. It is clear that the 
improvement is a function of the relative changes in the moment derivatives, 
iene end Mo . Both are reduced by the reduced flapping, and the mew rcor 
locations on the complex plane depend on which one is reduced "more’. In 
other words, at each forward velccity the suppression of lcngitudinal 
flapping will change the stability characteristics by scme amount, and 
that amount depends on the size of the flapping derivatives, 99/36 
and ee), = the relative amount of R tilt due to flapping, af 


and, finally, how much the rotor actually contributes to the stability 
derivative (e.g., the tail contributes the major portion of Hg at high 
speeds). These will all change with speed, and the plot of figure (12) 
shows that the relative changes in the derivatives for the "average’ 
helicopter are such as to lend credence to the flapp:ng feedback as speed 
is increased. Unfortunately, it has been showm that, even at the Nigh 
speeds where, from figure (12), the feedback seems most desirable, improve- 


ment in stability is only modest. 
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Collective Flapping (Coning) 


Near the beginning of Part II it was shovn that aa, and 40,have livtle, 
ie@eed Neelieibile, effect on the orientation of the rotor resulteank @oree 
weector, Rh, except wnen the helicopter is pitching. In the absence of 
pitching, the AQ, and “Ab, contributions to the H force are alnost 
equal and opposite in response to an angle of attack change. In a static 
sense it can then be said that coning angle has no direct influence on 
the stability derivatives. It is, however, acceptable to approximate 
EM PPoporvionality Of thrust and collective pitch, which in turn, in the 
absence of induced effects, is nearly proportional to the coninrg angle. 


It is approximated as 


WI Adz 


QO. im 





i = 


eucreait all of the lift comes from the thrust of the rotor, 





— — 


Zz Mere Se | 0 Oa, PL as 
a mM de Qsbreim Ou Qs dw 


It is now possible to write a feedback equation, using this simple 
@pproximation, that says the collective pitch of the blades is 


proportional to the coning angle 


A 6, Ale aa me: 


hg 





If the feedback loop is infinitely fast there is some value, Yes) 
Os 


of gain such that AG5= 0, and then oe = O which causes Zw = 0. 
Ow | 
In addition, now that the collective pitch can vary, there are changes 


in TPP tilt and R vector tilt due to the derivatives 


Oa and oo ots 
85 00: OGc 


fiitewitll be shown later, in the analysis of the actual feedback systems 
on the analog computer, that these changes in a, and a’ are of prime 
amportance when considering any kind of thrust feedback. They result 
in cross control (M@,) coupling and can significantly alter the angular 


accelerations on the helicopter. 


For the present discussion, only the effects of reduced Zv will be 
considered. When Zw = O the characteristic equation of the constant 


speed approximation becomes 


$(s-ta) -Ma = 0 


where the vertical force equation now reduces to the identity, AX=AQ. 
There is one unstable root for any positive Mex. Physically, this is 
the equation of motion of a weather vane. The sign and magnitude ofr 
ie depends on the position of the pivot point. If the pivot 38 gar 
enough back Mx is positive, and the weather vane diverges in rotation 


when disturbed, 





A reduction in Zw causes a decrease in dampirg and frequency of any 
aircraft. It also causes a reduction of initial normal acceleration. 

To the airplane, with its large negative My, a reduced dampirg and 
frequency of the short period oscillation is usually a small "cost" for 

an auto pilot designed to reduce Zw, and thus the initial lead factor. 

Per a helicopter with a positive My, however, a decreased Zw is rad onl; 
at the expense of an increased rate of pitching divergence. With a 
positive My it is mandatory to reduce the angle of attack quickly in order 
tO achieve pitching stability, but this means increased load factors. 

This is clearly one of those cases in which other restrictions must de 


considered when deciding which is best. 
Two Degrees Vs. Three Degrees of Freedom 


Before going on with the analysis of any particular feedback s:stem, 
it 2S of interest to see how the perfect flapping feedback systems 
considered above affect the long period oscillation, or tne whode of motion 
arising when forward speed is allowed to vary. Using the approximation 
derived in Part I, the roots for the 42000 lb helicopter neve been 
computed for both the unmodified and the modified cases and are shown in 
figure (13). The root locations show a modest improvement in stability. 
ine figure is not a conventional root loci, but rather the constant 
Speed roots were computed for the two cases individually, and a rcot loci 


flouted for increasing Mu. 


it vis Seen from the small change in long period roots Tmt There as 
very little effect of flappine in the long period oscillation, Prem vais 
observation it seems reasonable to consider onl; the consta:t speed 
approximate equations of motion for a simplified amalgsis. IU 1s proebaole 
that no improvement in helicopter stability characteristics caw be Made 
without improvement in the “AWzoOncde", and so this approximation will be 
used in the analytical portion of the remainder of this work. "Fae @weloc 


computer results will. shew the couplete dagrees of frealiel, ’iweve. 
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Oehmichen and Delta Three Piteh Flap Coupling 


Flapping feedback control can be implemented on a helicopter mechayically 
or electronically. ULinkages or sensors are used in some way to sense blade 
flapping and alter the pitch setting of the blades cyclically or collectively. 
There are two mechanical schemes which deserve special interest. The first 
is Oehmichen pitch flap coupling because, from an enal:ssis cf this mechanical 
device, insight into other more complex systems currently being prcposed for 
helicopter stabilization can be found. The secotfd is Delta Three pitch 
flap coupling because it has been considered as a possible means of achieving 
blade and helicopter stability for many years, but has not, to the author's 
Knowledge, been investigated in connection with the dynemics of the nelicopver 


as a whole. The two mechanical schemes are described below. 


fem Oehmichen. patch flap coupling the cyelic variation of piteheyaa, 
a blade "sees" (excluding any pilot inputs) is proporticual to the flapping 
angle of a blade (in the case of a four bladed rotor, but othervise this 
should be the angle of the TPP) 90° in front of it (in the direction of 


rotation). This is scen in the figure below for a four bleded rotor. 
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Positive Oehmichen coupling is defined here as: 
O. Ls By Coo W + Ais 3m YD =-H he A, Coos — by Su W 


as) 


in words, positive blade flapping resulting in a decrease in pitch of 
the blade 90° behind it (counter to the direction of rotation) is 


YL 


@etined as positive. The mechanical gearing ratio is the gainvH. 
The Delta Three hinge is the blade flapping hinge and when the 
hinge is cocked opposite to the direction of rotation it is said that 


o3= (I between the hinge line and the blade transverse axis). 


It is shown positive in the figure. 
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With positive Delta Three, wnen a blade flaps up it reduces its own pitch 


angle by tan $2. 


Thus AS = -TonS, OB 


Slit 


Qo - Ais Crap - Bis Ga = ~ten8 0-4 cry = bi rin W 


iets sce tet, with both schemes the collective pitch is decreased 
feewercviOnal to the conang angle, but that tne two devices differ in Ue 
mele component of feedback. The cyclic pitch due to the feediaick, is Alp 
fiemprepercional tO longitudinal flapping in the Oeumicnen system, but 


proportional to lateral flapping in the Delta Three sz stem. 


lnmereer  loOrecriuececrilvyely compare the two S¥svuems 10 15 CCBvesie7 aac 


Heparave the coning and cyclic components. This results in the feedosck 


equations: 
OBHMIC HEN DELTA THREE 
ADeo ee om 4 Ky Qo LSO-7 lo - J LO 
A Bis aad K AQ, ere ae ae Ab 
AAw = K Ab A Hig =e Keon 


where J = K = H, and J = K ¢ = tandz are the conventional, mechanical 


Ls 


arrangements Lor the two systems. 
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This separation in the feedback signals cay be accomplisted™oy suise@ole 
linkages or even by separate mechanical devices (reference 4), but, without 
regard for practicality, it is easier to think of the feedback es 
electrical signals originating from potentiometers mounted at the flapping 
hinges measuring blade flapping angles. In this way the signels may be 
used to drive servo-motors or hydraulic servos to vary the pitch of the 
blades at any azimuth position. The signals from different blades are 
easily combined and thus separation of the cyclic and collective components 
is possible. It is also possible witn this kind of arraigement to integrate 
the feedback signals and, with suitable sensors, to create feedoack sigrals 


proportional to the derivatives of blade flapping. 


Flapping feedback can be intreduced into the equations of motion on 
the analog computer very easily when the blade flappins variables are 
treated as they are in this work. On the analog computer diagram, figure (1), 
it is seen that the flapping variables are just different ccmbinations of the 
helicopter variables, each multiplied by a different constant. In actuality 
the different feedbacks are just unique combinations of angle of attack, 


airspeed, and pitching rate. As such, the feedback equaticns can be written 


as: 
Co 0= Osho = —) os Aw * Sot AU + Ode AQ 4 dus iS) + Ae N@ >» 
° ow Ou 0d OB: 
23 da, ~ : ox 


ee = el obi , wv Ob on 
ioe a K Nb, 2 | du Ow + Jy OY * SG NS) a aB abot Faas 





It should be noticed that the control variables appear on both sides 
Of the equations, ard there is therefore a closed loop effect which 
determines the actual "gain" of the system. This is illustrated in a 


simplified block diagram for the collective feedback. 
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iie@e porgion outlined by the dotted line can be replaced by @ sSimeie 
"block" which is the effective closed loop gain. One other olservatiol 
that should be made here is that the Delta Three and Ostmichen cyclic 
components are of opposite sign. Bis a, 18 increased wien the rotor flaps 
back (positive flapping) with positive Oehmichen, and “Ws Ap 1S decreased 
when the rotor tilts to the advancing side (positive flapping) with 


positive Delta Three. 
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Tne feedback equations can be rewritten in terms of effective 


e 
2 gains 


Which take into account the closed loop effects of the dependence of 


flapping on the control inputs. 


Sewn for the Delta Three system. 


Derivation of these effective gains 


Equations II (6) can be rewritten 


is 
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and substituting equations II (7) and II (8) into eacn other results in: 
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The Oehmichen feedback equations are analogous. Trey are: 
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Equations II (9) - II (12) are the complete feedback equations written 

an such a form that they can be introduced into the equations of motion, 
equation (1), as control inputs. They are multiplied by the respective 
control derivatives (e.g. MO, > Mp, , ete.) and transferred to the left hand 
side of equation (1). This procedure results in what ma; be thougit of as 
modified stability derivatives. It is shown for the constant speed memept 


equation modified by Oerhmichen feedback. From equation en 


-My aw + -Mé)d0 = My AW 4 Ma, AG. + Me, dBi 





which can be modified by replacing AD, and OB, dy equatiors TI (11) 
and II (12), resulting in: 


aM AwWa * (o> Me) NO = eas + Mg,|eqe T() | + Mg lew tor | 
Rearranging gives: 


/ / . / 
My Awa + (s-Ma ) Ae = fy elles 


where 
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The procedure is identical for all of the derivatives and for each 
proportional feedback system. It results in equations or motion of 
exactly the same form as without feedback, out with modified derivatives 
denoted by the (') prime. The writing of all of the modified derivatives 
is omitted here because little insight can be gained by looking at such 


complicated expressions whose terms may all be of the same size. 
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It is appropriate to analize each component of tke two feedback sjstems 
separately before considering the complete feedoack equations. There are 
three variations; lateral, longitudinal, or collective flapping may be 


used as the feedback "signal", independently of each other. 


Longitudinal (Oehmichen) Cyclic Feedback 


The feedback equation for this system is written as 


ABy = K 4a: A6,= O 


1) CASING ee CONN + Bou | 
Eee) \G Duy = AS Die “ 


which can be substituted into the equations of motion and then, regrouping, 


results in modified stability derivatives as follows. 


io 


7 ' [= 1 Oh 
My = Wien at Ma, Kk’ ee Zw ons eee >. 
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Oar Oa 
M ZW n : are always negative, 24 and are alwai\s 
Bi? > a d Oe a cL .LW cL | & 3 By ¢ dw c rd 2 my, 


positive. If Mw is positive, as for the wista@baliged helicopterg fan 


all of the primed derivatives are smwmller that the wiedified opes. Te 
resulting characteristic equation of mitra Is simiridy t mre Gi ae 
case previously considered, tnt with two aera BCR 





First, the "“autopilct” control is now being accomplished by tilting 
the control axis in response to a change in longitudinal flepping. Ir. 
the hypothetical @, = 0, it should be remembered, no provisicn is made 
for the tilt of the control axis - the flapping is "suppressed” in sdme 
unspecified manner, and any forces or moments dependert on flapping 
never even appear (see equation II (5)). In this actual Q,-> Bye feedback 
case, however, the control axis, ANF, is tilted as the “autopilot” respozse 
to flapping caused by a gust disturbance. The change of ANF angle of attack, 


BeCnde > Changes the H foree, thrust, and induces flapping. ATT of trese 


then change the moments on the helicopter. Tne,moment derivatives are 


changed by 


/ / 
IE (13) Moy =Mcy + Mak ay 


where Mg, can be expressed as: 
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By regrouping in equation II (14) to expose the implicit dependence on 
flapping, and substituting back into equation II (13), an expanded 


expression is found for a moment stability derivative. 
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It is seen that longitudinal (Oehmichen) cyclic feedback does not simvl,; 
relieve hub moments, but influences the moments on the helicopter due to 
each term in the above expression. It does not necessarily ccmpletely 


relieve hub moments, whereas the hypothetical a, = O case, earlier 


Coocicered, does. 
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The second point of importence is the occurrence of the cress control 
derivative, Zp,, in the expressions for the modified derivatives. This 
causes furtner differences between this actual feedback and the hypothetical 
ime O case, Ine cyclic control, By, is primarily a moment producing cGmezer, 
but vertical forces are also produced by this control input, and can cHarge 
Moe seharacleristices of the control response if the cross control derivataaes 


ZB, > is large. An expression for ZB, in initially level flight is 


a2 _ dT | a OAL eg 
—- = = — — oon q 
dB, wal ae | Od i 


but can be considerably simplified by neglecting all but the ay aa term, 


Viacom 42 good approximation in level flight. The derivative Chew geeomes 


| 
\l 


WA Big. ~4, = “UsLZw 


ZB, te neardily negligible in high speed flight where Zw 1S om the CGuG@ersor 


inasgy. When the approximation is substituted into the expression for Zw | 
7 my 
= 7 eG eee 


For all. positive K the derivative is reduced, and thus tie static Normal 
loads on the helicopter will be reduced, due to c§&clgie feedaack. Tue 
eross control causes little change in the "effective spéed"’, Ug', 73 


high speed flirt. 
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In figures (14) are shown the root locations for each data “elicopeer 
with (,-* Bis feedback gain of unity (top part of each figure), including 
the effects of Zp, . (The lines between roots are not conventional root 
loci.) It is seen that little difference between this and the A@,=O 
case shows up on the complex plane. Figures (15) and (16) shew analog 
computer responses for the 42000 lb helicopter with increasing gain, ¥, 
@Zp,= 0; and for varying Zg, @ constant K. Figure (16) snows resporises 


for K = .5, and K = 1.0, with 2g, at the value for the 42000 1b helicopter. 


The most important observation to be made concerning the analog 
computer responses is the change in sign and magnitude of the initial e 
response - Signifying a changed My . The helicopter pitches down 
dike an airplane with high cyclic feedback gain, and is thus staticall, 
Suge with angle of attack. It is also possible to see that the static 
normal load factor is reduced slightly wnen Zg, is increased from zero. 
Pijaily, 1c 1s apparent that cyclic feedback tends to stabilize tre 
helicopter in the longitudinal plane, but, for all of tke same reasons 
above, the same feedback causes a destabilizing lateral-lorgitudinal 


Coupling term. This is because the derivative, a or is positive, 


Causing positive lateral flapping in response to a vertical gust. “Weer 
positive lateral cyclic feedback is emploved, bAis a= KAb, , the Sisn 
of the feedback is such as to increase the longitudinal flapping irstead 


of reducing it, in response to a vertical gust. 


Lateral (Delta Three) Cyclic Feedback 


The cyclic component of Delta Three pitch flap coupling worzs tde 
same way as the longitudinal feedback just considered, but the sign ard 
magnitude of the feedback "signal" are different. The sign is different 
mecalse 2 rotor tilts with a positive b, due to a vertical gust, end 
positive Delta Three then causes a decreased longitudinal cyclic, A% sa). 
iitconissclearly the wrong sign for Longitudinal staticadtability. 1t 
feetne right sign, however, for blade static stability. Right fremp tee 
outset it is apparent that the positive Delta Three cyclic component will 
maerease helucopter longitudinal static instability, and snculd be avewaed 
from the point of view of this work. An investigation of negative Delta 
Mires 1s warranted from this observation, and positive Delta Three wall 


be disregarded until the collective component is included in the analjsis. 
ide procedure ©£Or analyzing this feedbaek as identical bo the 


longitudinal one. The modified stability derivatives tnat result can be 


written as (for the constant speed approximation): 


db, 
Mw > a Ke Ma, oa a = #4 We 23, a 


f ' dbi / / Oo 
Mé a ie - Ke Me, SO Us = Vg ~ Ks Za, 06 


When negalive Delta Three is used, all of tHe deriMmtives are Pétficed. 


Theiy magnitudes differ from the %Prictien fTeetlhbwck O8“LY YM LOS maT tudts 
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It 1S possible to compare the two tipes of feedback by comparing the 


flapping derivatives. It can be shown that 


Ody . qa lis | Dior 
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and 
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For normal values of ¥ (11-13) it is seen that 


Ole > On at = db 
Ou Ow O& Bye 


Paicweao 1Orward Speed increases, Ob; approaches Jai 

SS ow 
mie sitiiaraity of the two types of feedback is evidert, but the reductien 
in pitch damping is smaller, and the increase in static angle of attack 
stability is greater with the negative Delta Three feedback. This is seen 
in figure (17) where analog computer traces of helicopter motions with the 
different feedbacks is shown. The positive Delta Tnree arrangement is 
clearly destabilizing, and the negative Delta Three is clearly "better" 
than the Ochmichen feedback. In fact, the cyclic component of nepative 
Delta Three feedback alone would seem to be quite a good stabilization 


fevace ior a helicopter, on the basis of these results. 
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Collective (Coning) Feedback 


The collective feedback equation can be written 
NO, = — J Ade AGs =O 


and it should be noted that there is no distinction between positive 


Delta Three and Oehmichen because APis =O. 


When the feedback is introduced into the equations of motion, as vefore, 


and the approximation that 


Otte = DU» 
OW 08 


is made ( 200 arises only due to the change in local linear velocity at 
e 


the hub when pitching), the modified stability derivatives become 


~ 


( / So = 
Mee Muy ~ Mo, J ous Lie Exe ee ie 


Ma = He U. + U, 
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Hees see that, in addition to the anticipated result that @w' is 
meauced., Mw! is also reduced (made more statically stable) if 9, is 


positive. An expression for Ma. can be written as 


Oi OV ¢ fe Da. ate ebMs 12° Oa. 
005 A dd, t BOE 2 9, 


{ 
where — and oa, are both always positive. The sign of the 
= Q 


remaining term depends on the orientation of the R vector in trim. For 
the data helicopters the second two terms are always larger than the first 
(the first term is generally negative because a negative €results in a 
stable contribution to My), resulting in positive Mg. It can also be said 
that, unless the horizontal tail, fuselage, and hub contribute excessive 
nose up pitching moments in trim, Mg, must be positive and will become 
increasingly larger with forward speed. For the data helicopters Mg. 


increases by an order of magnitude from 60 to 150 knots. 


Gi 


Positive coning angle feedback thus reduces both Zw and tw’ for the 
helicopters represented in this study, resulting in a reduced maximum normal 
acceleration due to a gust disturbance, and an increase in angle of attack 
Beav2c Stability. These changes cause changes in the helicopter divigiigese: 
represented by shifts of the characteristic roots on the complex plane. 

As has been seen before, the reduced My will "improve", and the decreased 
Zn | Wawebevibtaacte the stability of the helicopter. The met cee) ae 


of course, depend on the relative sizes of M,, 4w, La. > and, especially, Mg. 





In the following figure this point is illustrated by showing the various 
yes aes | ! 
Mestt1oOns ©: the roots as Bw and Zw are reduced, and the resulting 


positions of the roots when Meg, is large and small. 
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The resulting root positions are thus very dependent on? 6, > the larger Mg, 


is, the greater the improverent in stability. 
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In figure (14) (bottom part of each entry) are shown the roct 
locations for the data helicopters when J =1. The very differer.t 
ways in which the feedoack has changed the roct locations for the 
four helicopters attests to the sensitivity of the coning feedbdacx 
to configuration. The 42000 1b helicopter exhibits the most pronounced 
change in that the helicopter is made stable with coning feedback. The 
analog computer responses of figure (18) show the dependerce of this 
helicopter on Mg. When My, is zero, the helicopter is divergent, but 
with Mg, set equal to its calculated value on the analog compute1, the 
helicopter response is considerably improved with coning angle feedvack. 
ies is8the only helicopter of the four, however, in which an appreciable 
improvement is shown. 

In order to determine the magnitude of the modified derivatives, 
J’ as a function of J, the mechanical gain must be known. According 


to the way this is defined, 


and thus at every different flight condition (different Qe ) 


2 


the actual "closed loop" gain is different. 





The function has the form (as do all of the "closed loop” gains in 


these feedback systems) 


i A 
[1+ %e | 


— 


and a sketch of this function is shown below 





B positwe 


B nec tive 


The function is linear thru zero, and if B is large, ah apprcoxjmetely 


linear region may be assumed for small A, and then A' = A. 


When B is small, positive A' is “gain limited", but negativ is 


miplitied rapidly toward infinity. 





Im’ this case, J = A and B= Ode , and it is thus Becessmar tc 
060 


investigate the nature of Oae/, 


From the expression given for 4g, 


= = = 


OM _ ¥ 4 en od ocr 
de; S) ea 3 ® 


In high speed flight Oy, approaches zero (negatively) due to the reduced 


mapereance OL induced effects. In the flight cohdition of interésu it ae 


8. | 


reasonable to assume 


Y 


| 


Le ae 


CY 


and a high speed estimate is 


“AC = 7 Con = 


From the sketch of the function for the "closed loop" gain it is secn 

that J' never gets-bigger than about (+) 2/3, but approaches negative 

infinity even with very small negative gain. When J‘ is, in fact, approacting 
negative infinity, Zw is approaching positive infinity and the roosts of the 


characteristic equation are approaching the zeroes of the equation 


(Sra)(Stb) - AZw (s- Me) 


when Mg, = 0. This would secm to prove that negative coning feedvacx is 
pood, but it is only because of an omission of the effects of the commm 
feedback on the effective angle of attack stability, My. nen Fe, is 
not exactly zero, Mw is approaching positive infinity and thus nigel éve 


coning feedback is impracticable. 





im “veppoense to 2 


vertical gust, the static response of the Mélicopter 


is 
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where, from equations II (9) thru II (12) 40s and AB 
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Figure (14) has already been referred to, but now will serve as a 
summary of the possibilities of longitudinal or collective flappirg 
feedback used individually. It shows that, for gains of unity, improve- 
ments in stability are not impressive, and certainly rot acceptable when 
the complexity of the mechanisms are considered. The analysis was not 
carried beyond gains of unity because it has been shcwn that there are 
ipade motion stability limits at, or near, unity for beth kinds of pitch 
flap coupling. References (3 and 6). These limits are the result of 
fixed shaft investigations of the dynamics "within" the rotor plane which 
have been assumed to be well damped and of high frequency in this study, 
and therefore negligible with respect to the motions of the relicopter. 
As the gains approach unity the results must be viewed with some 
apprehension as to their validity. An analysis whicr includes the blades 
as a degree of freedom coupled with the helicopter fuselage is needed to 
determine the exact stability limits. It is probable tnat these stavility 
limits will not exceed those found from the simpler fixed shaft study, 


and thus gains above unity will not be investigated. 


Comparison of the Complete Feedback Equations (Delta Three and Cehmichen) 


Combinations of the collective and cyclic flapping feedracks will new 
be analyzed with emphasis on comparing the effects of the lateral (Delta 
Three) and longitudinal (Oehmichen) cyclic components. The investication 
wet uwO parts; first, a static analytical analysis which revealauae® 
a particular feedback combination changes the heliccpter static stabdility; 
and second, analog computer responses showing the effects of the Treedbacks 


@n the Overall helicopter dynamics. 





The signs and magnitudes of the terms jin equations II (17) determine the 
size and magnitude of the left hand side of equations II (16) which can te 
thought of as modified static stability derivatives, Mr! and Zar, 
Since it is desired that static loads on a helicopter be reduced bz means 
of the feedback system being employed, and since Mo, and Zp, are positive, 
ele Mg, are negative, the required signs for equations II (17) are 

Yox%, (+) and SOY. (-). It is evident upon examination of 
equations II (17) that the two types of feedbacks will not be equally 


effective for corresponding mechanical gains. 


8> 
The collective "signals", 7S de , are different due to the 
e ; Blo ; O01 " T ; : 
mipherences between Vow and WwW, J and Jg. It has been shown 
at ow W at high speeds, and it can be shown from the 
definitions of J" and Js" that J" will always be larger than Js" for the 


Saieme. In MOSt cases these two values will be relatively close, HGweéver, 


Sioeemeemch gain, J, for K = -1, 0, +1. The negative valuesme? SO/Ang 
correspond to reducing static normal loads and improving angle of attack 
stability (Mg, Oia ). From the figure it is seen that the cyclic 
component of feedback does not influence the collective feedvack signal 
poprecigbly. This is because, in equations II Glas ke 09% is the 
fietianeny term, even for large K or Ks. Among the four data helicopiers, 
however, the 42000 lb one exhibits a marked difference in magnitude between 
J” and J," at high gains due to a particular configuration of small Oda of 
(large K') and large ee F 56, resulting in a strong amplification 

of J" within the range of gains being considered (-)<3<+)). “his is 


shown below. 
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J never becomes infinite, but the sensitivity of Ag (Oehmichen) 


Sra OA UNIT} rye oo Gare. 


Momence value Of gain, J, increases markedly as J increases. In Tactjmaau: 
this helicopter high gain responses could not even be ovserved on the 
analog computer due to this high gain sensitivity. A plot of the static 
feedback signal, AQ /usg vs J (corresponding to figure (190), for this 


"sensitive' helicopter is like the sketch below. 
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Only near the limits of useable gain, and only for this special 
configuration does this phenomenon occur, and so, in general, the static 
collective feedback signal, and thus the initial normal acceleration due 
GO a gust, is relatively independent of the cyclic feedback open loop 


gain. 


In order to evaluate the effectiveness of the static cyclic feedback 
Signal in changing the helicopter static stability, it is helpful to 
expand all of the terms into the more fundamental rotor aerodynamic 
derivatives. Before this is done, however, it should be noted that the 
equations II (17a) are identical in form, of differing sigzi, and thet 
when J = O and K and Kg are positive, the Oehmichen feedback signal is 
Stabilizing and the Delta Three is destabilizing. As J increases, Mowever, 
tiewoatenals decrease and, for some J, change sign. After uch algebraggac 


expanded expressions can be written as: 


OKHMICHEN 
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The static cyclic feedback signals will thus change sign when 





OEHMICHEN DELTA TH#HS 
VA 
te oe i (+38 | eee 
f+ 2%) e fe £24 


In the absence of induced effects, as a high speed approximation, and 


for %=12 the feedbacks change sign when 


OHMIC HEN DELTA TEREG 


y/s 7 RR = 


It was shown before that J' is gain limited at about 2/3, but for J<1 (bdlade 
flapping stability limit), J' never will exceed about 1/2, even at the 
hiehneaw speeds. With this information, it 1S predictable that tie Bele 
Three eyclic feedback component will always be statically destepilizing, 

and that at high gains and high forward velocities the Oekmicher c7clic 
feedback component will change from stabilizing to destadilizing. This 
prediction is fairly independent of anything but airspeed. As the air- 

speed is reduced the induced effects boost the destabilizine limit of J', and, 
Giethe same time, reduce the actual value of J’ so that alt low Speed Gat ee- 


memchanee In sien of the feedbdack. 


The prediction is verified in figure (19a) where the statie eyelic 
Feedback signal is plotted for the Oenmichen and Delta tee feedbacxs. 
it is seen that, for the conventional mechanical arrangement of Ocm—nctie 
Fiteh flap coupling, the static feedback signal changes to the ‘wromge' siggy 
as the mechanical gain is increased to near unity on tre 24000 1b helicopter. 


The same result is found for all four data helicopters, te crossdvet pods 


being close to that shown in fcume (19), but the slepe of the clirwe vars 


with configuration. The sane cenfigmuration (42000 1 Nelicimites’) ti) shtarei 
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The reason for this is the increased importance, in the cyclic feedback, 
of the increasingly larger collective feedback signals, as the mechanical 


Bain of the system is increased. 


As a summary of the findings of this static analvsis, it can be simply 
said that; 1) Positive collective feedback reduces static loads (r.ormal as 
well as angular) almost independently of the sign or magnitude of the cyclic 
feedback, and Mg, plays an important part in how "good" this t;pe of feedback 
is. 

2) Positive Delta Three cyclic component always increases the static angle 

Of avbvack instability of the helicopter, and negative Delta Three is alveve 
stabilizing. 

3) Oehmichen cyclic component will improve static stabilit; for most of the 
range of useable gain, J, but, at some high gain and high forward velocity, 


mall begin to vitiate the static stability. 





PeeCpmriese Sualic results arises a very interesting helicopter 
stabilization technique - negative Delta Three hinges on the mein 
fmouor. it has been shown that this technique improves static stamp, 
but two questions must be arswered before it can be seriously corsidered. 
First is the question of blade motion stability (this work cannot uncover 
blade instability because of the assumptions made avout the rotor plane) 
which immediately comes up because of the nature of the piter-flap coupling 
(the pitch of a blade is increased when it flaps up) is intrinsically unstable. 
This matter has been investigated in reference (6) where it is shown that 
"The use of negative Delta Three introduces the possibility of flapping 
divergence but this is not a problem for normally required values of 


negative Delta Three’. 


The second question concerns the increase in Zr' associated with 
negative Delta Three hinges. According to the analysis of this worn, 
in which only lowest order rotor plane dynamics are considered, ore 
fa7ecuwor negative Delta Three is to increase static normal pust Loads: 
figs is undesirable, but the validity of neglecting an; higher order 
dynamics (gust alleviation as defined in Appendix IV) is questiorable, 
especially with feedbacks such as these, and it is douodtful if the adverse 
effects of negative Delta Three are as pronounced as they seem in this 
analysis. A more complete analysis, however, is needed to verify the 


above assertion. 


The dynamic responses of one of the data helicopters with various 
feedback component combinations are presented as figures (20 - 23). 
As expected, those feedback combinations which improve static stability 
(My ) are also the ones which improve the dynamic gust responses. As 
such, negative Delta Three (figures 21 and 23b) feedback is cleerly the 


most promising feedback of tne various ones investigated. 
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Figures (20 and 21) show that gust responses are improved for: 


OEHMIC HEN DELTA TREEE 
K = constant, J increasing J = constant, X decreasing 
J = constant,(K = some intermediate (increasing negatively) 


value before crossover to negative feedback) 


The analog computer responses are revealing and, in the case of the 
negative Delta Three feedbacks, the improvements in gust response are 


seen to be excellent. 


All of the results above are in agreement with, and may expleir, tke 
surprisingly successful flignt test results briefly reported in reference (6). 
In that test a helicopter was fitted with negative Delta Three hinges 
(8,= -39°, Ks = .81) and flown to its power limit. The investigation was 
directed at rotor blade flapping stability and amplitude,and so no data on 
iMemwetootlitgy Of the helicopter is presented 1% “hat report. 


Oehmichen and Delta Three Variations; Rate and Attitude 


Several helicopter feedback concepts tnal have been introduced in the 
literature can be quite adequately described by variations of the Cehnic’ten 
and Delta Three feedback components when only the macroscopic dynamics of 
the helicopter are being investigated. In reference (9) is proposed a thrust 
and tilting moment feedback system, in reference (13) is preserted the Lock- 
Mecareyro system, and in reference (4) are presented various stabilization 


schemes. All of the above fit into this category. 
The proposed thrust feedback can be pretty well descrived b} the 


integral of the collective component of Oehmichen feedvack, hO,~ -rf bas : 


where it is an adequate approximation that 4 ey Cs Ve 


8? 








From experience with the non-integratéd feedback (460, = -J Aad.) it is 
predictable that M@, plays an important part in tne effectiveness of the 
integrated feedback. As can be seen from figure (24), when Mg, LSegerS 
this feedback scheme vitiates the helicopter response, but when Mg, i 
positive, increasing the gain increases the frequency ard decreases the 
time to double amplitude of the response. This is similar to the result 
found with the non-integrated collective feedback where the static atgle 


of attack instability was decreased due to the cross control derivative. 


The integral coning angle feedback (thrust feedback) dces not reduce 
static normal. loads, and from figure (24) evidently reduces the effective 
Zw of the helicopter out of phase with tne angle of attack variation 
(the integration makes the loop too "slow"), and thus only improves 
helicopter dynamic response if the cross control derivative, M6, 


is large and positive. 


Figure (25a) shows the results of integrating tlhe cyclic cormponelt of 
Oehmichen feedback, Bys = x fod, dt as proposed in reference (4). No 
Static stability derivatives are modified with this sciteme ad teLacopter 
dynamic gust responses are not improved. Again, the intefratiew mates the 


feedback too "slow". 


Combining the two above, or integrating the entire Oehnmichen feeddack 
signal, results in the gust responses of figure (5b). As with tre 
individual component feedback signals, the loop is too "slay" amd dees rot 
improve the response of the helicopter. This particular combination is 

- ‘ 


the same as that proposed in reference (9), the Thrust aiid Tilting “owetit 


feedback. 


The one variation of the Oehmichen feedieack that deserves special 
attention at this time is the Lockheed gyro systeth. This s¥stem Was eth 


devised in order to relieve dyn@iiic rotor loads @id to sUBAalfve ie Nelé- 


copter, as a whole, inf the presevce of gusts. The fdrel fears cagies Fe 
treated within the seoyt of this work, ad co the WHS!) vsed Per wre seer 
herein id@ilives the VWerhines of tbs ety, ame oO) Alesse neDes 4 


second a] gamed femtaime 





In the ideal system L = 0,Y=1 (in the nomenclature of referer.ce 13) 
and the feedback reduces to the integral cf Oehmicner. plus tne intezral 
of pitching rate, for the purposes of this work. A block diagrew cf tre 


idealized system is shown below. 
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It 1s not possible to perform the seme static feedhack amwalwsis as 
was done on the Oehmichen and Delta Three systems because of the integration: 
which does not “let through" any static feedback. A feedback equatlo car. 
be written, however, which resembles the ones of the previous sections. 


From the diagram: 
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"sweep constant” = F(P*-1) in reference (13) 
& 


where K 4 


linkage between gyro and swash 


and Co 


Except for the added pitching rate tezm, and the lag time constant 
associated with the feedback (s in the denominator), this Gyro feedback 
equation resembles the Oehmichen equation. Typical values for Ky ame 


Q- 1/2 and Co = 2/3. 
cae” 
It was found before that Oa = Oey. rt (3 ee 
I U i | 5 x AA so | 


so that even with large Ky there is a long lag time, or the feedback 
loop is "slow". Because of this long lag time the gyro feedback dces not 
aporeciably alter the motions before tre angle cf attack convergence, sd geetce 
Waemiormna! acceleration response, is completed. Since there 1s "6 Coli ga 
feedback there is only a small alteration in the norwlal force equataon 

(28,88, ). The lagged rate feedback acts as a torsicnal spring, but 

Sloe the rotational motions only slightly. For all cf tmese reacam 

the analog computer responses (figure 25c) for the idealized Leckheed 

gyro feedback system are sorely disappointing. The responses are for 


Cy = 2/3 and increasing Khe 


As a standard for comparison of all of these feedback lechmigies 7: me 
have oeen looked at, it is now appropriate to show Lhe Bust respamse of a 
helicopter equipped with conventional SAS. Figure (26) presents the analog 
computer responses of two of the data helicopters modified by rate dima 
attitude feedback - the last two strips of the figure are for values of 
teat 5: and R with which the production models of those helicopters are 


presently equipped. 


The holicovter responses are, at first glawce, quite "goml’. [1 st@ns 
that there should be little reason for a pyloy lo c@igmlaey.. Recipe iy da. 
however, it Bhcones nechRkae Lo 1@e Wer. cramely ats PeMreltwesS ye al 


The modal weed for chlai blew UO 1. ie ia < PaLL! % Le 





on here, as this work is not designed to provide very accurate detailed 


models. 


In the figure it is seen that the acceleration responses, weile meawaly 
damped, exhibit considerable static (time = 0 a values @nd over Siigec. 
Because this feedback does not alter the static response (2 ), gust induced 
loads (normal and torque about the center of gravit:) are unchenged. This 
means that the pilot will still feel the same local "eg" (i.e. the normal 
load, Zw wg, plus the upward angular acceleration times his distarce from 
the C.G.) response as he did without the SAS. Even though the motiors are 
then dissipated quickly, the pilot is still subject to the same maximum 


forces on his body with this feedback. 


Tt is evident, from the point of view of this study, that tie 
improvement in angular static stability is the prime determinant of the 
quality of a feedback technique. As such, no flapping feedback system 
with an integrator in the loop, and no body mounted SAS which does not 
sense angular acceleration will appreciably improve the helicopter pilot's 


meron Of the flying qualities of his craft. 
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CONCLUSTONS 


Parts I and II, being somewnat separate from each other in material 
presented, warrant separate conclusions. A summary of the fiddings i” 
Part I is given on pages (28 and 29 ) as the physical interpretatio:. of 


tiewcliuibarities and differences found in that part. In addition: 


1) Airplane and helicopter Zw's are found to be quite nearly equal 
when comparing aircraft of similar high speed (150 kts) design corditions, 
This leads to the discovery that the normal acceleration gust response 


of the two types of aircraft is similar. 


2) The size of' Mg and the sign and magnitude of My are different 
Persuie UuvO types Of aircraft, the helicopter being on the unsteblescide 
im both cases. This leads to the realization that the rotatiotal Motions 


are of prime importance in helicopter gust response. 


3) The influence of the fuselage, tail, and rotor flapping Hinge 
offset (or non-articulated blades) becomes more importunt than, and in many 
cases overshadows, conventional” heliccpter stability characteristics 
in high speed flight with the "typical" helicopter characterized hevei’.. 
This is to say that rotor flapping beccmes less important in the presenesc 


Of these Other factors. 


kh) Of the various stabilization techniques, involving rotor flapping es 
the basis for the feedback signal studied in this work, negative Delta P¥free 
(positive pitch-flap coupling) shoved the best (and only acceptable) resylts. 
The static stability is improved consideravly with this tfpe of Temameck. 
It is considered that this finding warrants further detdbled stud. i” Ge 


ered. 


5) Integrated feedbacks, the sigruls for whieh ordinate as rater 
flapping variables, ure too "slow’ to itiprove the gust rvs che eaierisstics 
of the helicepter. Petause cf lhe abd gates 1: 1 AOD, Ure fm eo 
in stutlic stelility oom CORE) se jae ie Re AAO APR 


response. 
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FIGURE (5). Three degrees of freedom analog responses 
for basic helicopters. a) 42000 1b, 1) 35000 lb, e) 24000 lb, 
ad) 18000 1b. 
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Longitudinal (Oehrichen) flapping feedback, 
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FIGURE (16): Longitudinal (Oehmichen) flapping feedbeck , 
Bs @ K a,, for a) K = 0.5, varying cross control derivative, 


Zg,* -05, -25, ») 2, ™ value of 42000 1b helicopter, varying 


K + 4 6.23 20. 
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FIGURE(17): Cyclic component of Delta Three feedback for 
configuration similar to, but not exactly, the 42000 1b heli- 
copter, a) positive (conventional) , Kg = 0.1, b) negative, 


K¢ = “1.2, -~.15, 1.0. 
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(a) 


Coning angle (collective flapping) feedback, 


FIGURE (15): 


b) Mg, ™ Value of 


a) Mg,2 0, J # .25, 0.5, 1.0. 


6, = -J a, 


42000 lb helicopter, J = 0.5, 1.0. 
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FIGURE (20). Oehmichen pitch flap coupling on 4200 





helicopter, a) Effect of coning component gain, J = 0, 0.5, 


1.0, @ K==0.5. b) Effect of cyclic component gain, K = .25, 


79, @ 3 @ 0.5. 



































uffect of 


Delta Three pitch flap coupling, 


FIGURE (21): 


sign and magnitude of cyclic component on configuration 


Ke = *.25, 


similar to, but not exectly, the 42000 1b helicopter, 


25, ~1.0, @J#0.5. 
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FIGURE (22): Oehmichen pitch flap coupling on 42000 1b 


helicopter; K # J # .25, 0.5, .75. 
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FIGURE (23): Delta Three pitch flap coupling for configuration 


similar to, but not exactly, the 42000 lb helicopter, a) Positive 
(conventional) J = Ki = .25, 0-5. b) Negative J = Ky, = -0.5, 
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FIGURE (24): Integrated coning angle feedbeck, @, = -J J a. dat, 
ao) Me, ® value of 42000 1b helicopter, J = .25, .15, 1.0, 


bv) My = 0, a .25. 
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(a) 
FIGURE (25): On the 42000 1b helicopter... a) Integrated 
cyclic component of Oehmichen feedback, By. = K{ a,; K = .05, 


5, -75. b) Integrated Oehmichen Feedback , Bic ® Kf a, dt, 


2 


Oo =- T{ng dt, J=h = .25, 0.§, c) Idealized Lockheed 


Gyro, C, = 2/3, hk, = 0.1, O.4, 0.6. 














VALUES ON ACTUAC #ELICOP7ER 


(a) (b) 


FIGURE (26): Rate and Attitude feedback; (+t,81)6 B,, = R(t, s#1) 5 0, 
a) 42000 lb helicopter, R = .261; 1)%,= 0, ti™ .26, 2)1, #0, 
t © 2.6, 3) 0.22, = .26, b)t, = .22, t,= 2.6, bd) 18000 1b 


helicopter, R= .55,% 3 .1, %* .37. 
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HELICOPTER DATA _ 


The data contained in this appendix is taken from unpublished data on Gdarce 
Waren nieh performance, single rotor helicopters. Forward and ait cencer so, 
pravity positions, heavy and light gross weights, and three different configuratichs 
pieerepresented. All of the helicopters have certain characterisvics ice 
and these characteristics are taken as representative cf any heavy, high speed, 
single rotor, helicopter. They include large flapping hinge offset (“/R = 5%), 

5 and 6 bladed main rotors, high disk loading, and relatively large harizontal 


teaiis, 


pome Of the data did not cover the flight speed of interest iI'this we. 
SO it waS interpolated to extend 10 knots and 15 knots in two cases. The 


following table shows some of the pertinent information ccncerning the heli- 


CopLers. 

GROSS CEG, MAX SPEED TATI e a 
WEIGHT LOCATION OF DATA (KTS) SIZE(s; 1, ) R A 
42000 FWD. VO 1040 >. on 9.3 
35000 AFT. 170 1670 5.58 8.2 
21,000 wD. UG 1670 5.58 a2 
18000 AFT. eS O55 3.2% ones 


The data was reduced from a spuce fixed axis to a stability axis «rd 


| dimensionalized. At each speed the high, low and average of each st#ui lit 


derivative was found and the informtetion is prese:r.ted in the fignpess “Te 
lines thru the data points represert "Lepiedl", or "avermmé” stat et! 
derivatives, and the tremam@ of virigtig¢ will Corpwabd Syeen. 
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APPENDIX ITI 


AXIS SYSTEMS AND AXTS COUNVERSTON 


There are three types of axis systems commonly used in aircraft 
stability analyses, each having its own peculiar advantages. The 
Conversion of data Irom one system to another is usually required when 
data comes from different sources or wnen the analyst is mcre familiar 
With one system. The axis systems are named by the orier.tation of the 
Peeeose ss Harst, there is the horizontal, or Garth fixed axis system. 

In this system the relative wind velocities are resolved into horizontal 
@idevertical components. Second, there is a wind axis system in which 
the X axis is directed into the relative wind at all times. In both 
of the axis systems above, an aircraft attitude change (e.g. angle between 
wing chord and horizon) results in a new orientation of the X axis with 
respect to the aircraft. This changes the velocity coniponents with 
respect to the aircraft, and moments and forces result. This results 
in non-zero stability derivatives like 

ax , NM, 22 

vax) oD “0 
It is important to realize that these derivatives arise only because of 
the choice of axis system. Aerodynamic forees do not depend on pitch 
altitude of the aircraft body (e.g. wing chord) with respect to the 


horizon. 


The one axis system left is the body axis system. As the name 
implies the X axis is fixed with respect to the bod; (e.g. a wing chord 
or a rotor shaft) and al] velocities are resolved parallel ai.d 
perpendicular to the body. With such a system an attitude chénmee does 
not change the velocity compotients with respect to the bGa@y" aud Merce no 


aerodynamic forces arise. This system is cheracltergaged Wye the alseoigee 


oui ok Garner HE der i Vai ves 1 i GQ : but rebuld i 3) } ral bi { O a a. 7 1s * - Ww ™ 
trying to interpret ppeeicall]. tyeeTorccs*tm 16 Wie 
All? J 








This is because the rotating (in the case of a pitching or rolling rate) 
Peeremeve System introduces inertial, or centrifugal, forces in order 

to satisfy Newton's laws. Care must be taken to remember thet these 
inertial forces arise only because the interpreter is in a fixed coordinete 
system and the aircraft in a rotating system. As such, the centrifugal 
forces are not really forces on the aircraft at all, but rather a may, fcr 
the interpreter to satisfy himself that F = ma without having to switch 


to the rotating axes of the aircraft. 


A special case of the body axis system is when the X axis is initiall); 
aligned with the relative wind, and then fixed to the body - all at time 
zero - but then stays fixed with that orientation in the body during air- 
craft motions. This further simplifies the equations of motion by 


eliminating a wo term in the Z force equation. 


CONVERSION FROM SPACK FIXED TO STABILITY AXES 
In one way of converting data from one axis System to another the 
velocity components in the first are resolved into the second and sub- 
stituted into the equations of motion, expressed in the variables of the 
first. Regrouping results in equations of the same form but with different 
coefficients. These new coefficients are renamed the stability derivatives 
in the second axis system. This is demonstrated by the following exaple 


of converting from space fixed to stability axes. 






= es _ es 
X space fixed 


ZL SAACE fraed 


Ail ¢ 





From the diagram of the velocity vector eccmponents comes the 
relationship 


8 


= Uy We 5 


W= Wz t+ uy,® 


where u, w are body axis velocities and u,, wz 


are horizontal (space 
fixed) velocities. Expanding, for small perturbations, gives 


Au = Aux — (Wz, AO + @,dwz) 
mw (1) 
AW = SWz + (Ux, 450 + OB AUx) 


li@ meee ly axes are déeiined by the initial angle of ettacze beim 
ewe SO, 1Or Stability axes, 


Oe SG ERS Ne 


i imeoaddition, initially straight level flight is comidered vien 


The velocity vector resolution, equation A (1), then becomes 
AW = AUx 

A (2) 
Aw ee DW x Ux, AG ae AW 2 + 1h Ue 


eer 7etorce equation may be written, in a space waxis. as 


elite Oe EY oF 2 AU = oF. SWz 4°? - 
Oux OUx 


=O 


Pulelal 


, 
) 





Substituting the equations A (2) into the above, 


OFe nu, + Mkean + Wray 4 We pw 
why OU x OW ows 








4 dre _ dbz ert abe] . Fz ‘SQ =06 
06 dkz d6lu, Xa Ju, 


The first bracket term reduces to 


SING) = OPE AS 
oY 


which is the centrifugal force that arises due to the chafize in axis 
systems (from stationary to rotating). The second bracketed term is 
identically zero (this is emphasized by the indication of what variables 
are held constant when those partial derivatives were teken). The cther 


coefficients are now renamed in the body axis as (after dividing them by m) 


Mor 2asu + Zu du + ean + Gene = Osorno 


The terms Z, and Z; are due to unsteady aerodynamics and are zero im tis 


work, 


Alii} 


Another method that can be used to convert frcm one axis to 
another may be more straight forward and can provide more insight irto 
the nature of the two axis systems. From the force balance diagren 
below come the expressions for body forces in terms of horizontal fcrces. 


(The moments are equel in both systems) 


= Fy Ca F26 
ji = be + ee 
M= IM 





Each body axis derivative is then found by applying the cNain rule as is 


demonstrated for Xu below. 


dX a dV x OMx a Obs aux (Q,) = ee dO_ OWUx 
OUlwe Uy OU Ux OW ee Oy oe 


, OFR due - ke O12 (g,) - Fy 30 DUr 


OW 2 ou OAyz dU Oke OU 


Alli § 








From equation A (1) comes 


eae OUx = © 
du OW 
OUz — Ne OU — - ©, 
OU OU 
And in the equation A (5) 
9 28 ay» 0 = Het Ke 
OU ~ OWz 
Rewriting equation A (5) for any bod; axes 
K 2 we. Wey, + yf ode - ee — 
© > tae 
OU OU x Oils U2 OUx OU x 
fanelly, tor stability axes 
ae ©): 
and 
JX = Okx -- Xo re 
Ou OU OU2 
re) Ot OWUn 


AGL & 





The same procedure is carried through for every derivative in the bed;, 


axis equations of motion. The result of all the conversion is: 


STAB ITT yeas SPACE AXIS 
7 |p a Ee ¥o[Fxuz eg Bo Fzuz | 
Kaos FxU 2 - Pol hue. ie ¥oFzu x ] 
4. = Fu oni %o [Fxu a Pau 27 Oo Fu Z ] 
Zw Fzu 2 + ‘ [Fone bean z + By hou va 
Mu = Muy - % My 
My = Mwa + Yo Mu o 


MA = M6 


pa We 





AreaypiA iV 
GUST ALLEVIATION 


for the Sawple cirerat, of this wore the initial (at time = 0 ry load 
factor is instantaneous and proportional to the gust magnitude and Zw. In 
powder Uieuhcaleielicopter Or airplane, this is not the observed 
response. The differences between The actual and the model responses are 
attributed to and defined herein as gust alleviaticn. Tne definition used 


here is in contrast to another definition that would interpret this whole 


work as a study of means by which to alleviate the magnitude of a relicopter's 


response to gusts, and thus considers artificial stabilization as gust 
alleviation. Many parameters may influence the mesnitude and dynamics of 
gust alleviation. In reference (8), for example, is shown the derivation 
fer a peust alleviation factor, or a corrective factor to mullipiy Cages 
Bie stau1c value oF 77, Which should be used to Compute the Initiale 
factor for a static model (e.g. rotor bolted down in a wind tunel) due to 
a more complicated inclusion of the rotor dcwnwash. In reference (5) is 
shown the result of a detailed computer program. used to compute gust 
alleviation factors for the initial and the dynamic response by ircludinge 
unsteady aerodynamics and aero-elasticity to compute an effective value 

of Zw that varies with time and the motion of the helicopter. <A very 
detailed computer analysis such as that used in reference (5), hovever, 
may or may not have included every cause of gust alleviaticn, and it is 
hard to gather any insight as to which are the important causes end their 
davect vetiects. ~Tne only conelusion that Camm oomiy be drag 1S° fiat eter 
helicopters encountering sharp edged gusts, the elfective value of “4p is 


dynamic and may vary from the static value. 


In the figure below are shown examples of what load factor respotises 
could look like for the static, zero degrees of freedom medel if gist 
alleviation were considered. (Note: these are not computed or predicted 


responses but are drawn for illustration only). 
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The value at time = 0 is shom as Zw for the linear model, aildnhz 
for the kind of static analysis presented in reference (8). The dotted 
line is descriptive of the results that would be obtained by extending the 
static analysis in an iterative manner - computing the thrust at each 
increment in time, finding a corresponding downwash velocity and using 
that to compute the thrust at the next point, and so on. The remaining 
dine illustrates the kind of dynamic response that mignt be found to 
ensue with an analysis such as in reference (5). In accord with 
assumptions concerning a quasi-static analysis, and neglecting unsteady 
aerodynamics and rotor plane dynamics in this work, a constant linear 
value of Zw is used corresponding to the value to the right of tne 


hatched area on the figure. 


In the static model the only aerodynamic variable that the aircraft 


responds to is AWa because MO =BW, =O. 


As degrees of freedom are added the motions of the aircraft couple 
with the gust velocities to cause aerodynamic variables that are functors 
of time and the aircraft motions. As the frequency of the inctioNs increase 
from zero (as in the static case) to the full three degrees of freedvr. 
motions, the errors of neglecting alleviation increase ad Jecmme obscured. 


The same is true if the frequency of the pust distumWalces is ecygemen. 
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It is helpful to think of gust alleviation as a dynarsic transfer 


function to replace the Zw and Mw potenticmeters on an analog computer. 
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At very low frequency the transfer function reduces to Zw, as the 


diagram shows, but at a higher frequency there may be a dynanic re: SiS. 
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Two ways of conducting an analysis are contrasted here: A detailed 
analysis in which every conceivable cause and effect is included for 
analysis, and a simple analysis in which only the most important, or 
tipsl Order, Ciie@eus are considered. [nis study is of the lavuter 
elassification. Until the first order, important causes of alleviation 
can be determined and modeled simply, however, there is little justification 
for the inclusion of any alleviation in the simple study. Other irvestigations, 
including flight testing, must be performed to determine which are the first 
order causes and what are their important effects. The secpe of this paper 


does not include such investigation. 


A pertinent study would be to compare the results of simple analyses 
such as this with full scale flight test gust response data. This could 
lead to some insight as to how alleviation really influences the dynamics 
as a whole. It might then be instructive to construct different filters 
on the analog computer and put them between the gust inputs and the 
summing amplifiers to try to reproduce the flight test data. The 
characteristics of a particular filter that did the job might give a clue 


as to what the real causes are. 
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